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Abstract — Consider the task of sampling and reconstructing a 
bandlimited spatial field in using moving sensors that take 
measurements along their path. It is inexpensive to increase the 
sampling rate along the paths of the sensors but more expensive 
to increase the total distance traveled by the sensors per unit 
area, which we call the path density. In this paper we introduce 
the problem of designing sensor trajectories that are minimal 
in path density subject to the condition that the measurements 
of the field on these trajectories admit perfect reconstruction of 
bandlimited fields. We study various possible designs of sampling 
trajectories. Generalizing some ideas from the classical theory of 
sampling on lattices, we obtain necessary and sufficient conditions 
on the trajectories for perfect reconstruction. We show that a 
single set of equispaced parallel lines has the lowest path density 
from certain restricted classes of trajectories that admit perfect 
reconstruction. 

We then generalize some of our results to higher dimensions. 
We first obtain results on designing sampling trajectories in 
higher dimensional fields. Further, interpreting trajectories as 
1 -dimensional manifolds, we extend some of our ideas to higher 
dimensional sampling manifolds. We formulate the problem of 
designing K-dimensional sampling manifolds for (j-dimensional 
spatial fields that are minimal in manifold density, a natural 
generalization of the path density. We show that our results 
on sampling trajectories for fields in can be generalized to 
analogous results on d 1 -dimensional sampling manifolds for 
d-dimensional spatial fields. 



L Introduction 

A. Problem description 

Consider the problem of sampling a d-dimensional time- 
invariant spatial field f{r) : r G M'', where r represents a 
d-dimensional spatial location. If / is bandlimited, results 
from classical sampling theory (see, e.g., ||3|, ID) provide 
schemes for sampling and reconstructing the field based on 
measurements of the field at a countable number of spatial 
locations, e.g., points on a lattice or a non-uniform collection 
of points like the one depicted in Figure [T(a)] The performance 
metric used in designing such sampling schemes is the sam- 
pling density - i.e., the number of sampling locations per unit 
spatial volume. Such a metric is motivated by the fact that 
one typically employs static sensors to measure the field at 
their locations and hence the sampling density is equal to the 
spatial density of the sensor deployment. 

Some of the results presented here were pubhshed in an abridged form in 
(1] and d. 




(a) Static samphng on points (b) Mobile sampling on a curve 

Fig. 1. Two approaches for sampling a field in 



The scenario is different in some practical cases. Consider 
for instance the problem of sampling a d-dimensional spatial 
field (where d = 2 or 3) using a mobile sensor that moves 
along a continuous path through space and takes measurements 
along its path. An example of the path of a sensor moving 
in is shown in Figure |l(b)| In such cases it is often 
inexpensive to increase the spatial sampling rate along the 
sensor's path. Hence it is reasonable to assume that the sensor 
can record the field values at an arbitrarily high but finite 
resolution on its path. The objective now is to reconstruct the 
d-dimensional field using only the values of the field at closely 
spaced points on the path of the sensor through R**. For such 
a sampling scheme, the density of the sampling points in R'' 
is no longer a relevant performance metric. Instead, a more 
relevant metric is the average distance that needs to be traveled 
by the sensor per unit spatial volume (or area, for d ~ 2). We 
call this metric the path density. Such a metric is relevant 
in applications like environmental monitoring using moving 
sensors 13], 161, where the path density directly measures 
the distance moved by the sensor per unit area. This metric 
is also useful in designing fc-space trajectories for Magnetic 
Resonance Imaging (MRl) [7], where the path density captures 
the total length of the trajectories per unit area in fc-space 
which can be used as a proxy for the total scanning time per 
unit area in fc-space. 

Now consider a different related problem. Suppose we 
want to reconstruct a 3-dimensional bandlimited spatial field 
using measurements of the field along 2-dimensional surfaces. 
Such a scheme is employed in applications like Transmission 
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(a) Sampling lattice (b) Sampling lines (c) Sampling planes 

Fig. 2. Three approaches for sampling a field in M.^: sampling on a lattice, 
along hnes, and along planes. 



Electron Microscopy (TEM) JU and MRI 19]. In these cases 
it may be relatively inexpensive to increase the sampling 
resolution on the 2-dimensional measurement surfaces but 
it may be more expensive to increase the total area of the 
measurement surfaces. Thus it may be reasonable to assume 
that the measurements reveal the value of the field at arbitrarily 
high resolutions on the surfaces. The objective now is to 
design measurement surfaces that admit perfect reconstruction 
of bandlimited fields and are simultaneously minimal in area. 

Motivated by such problems we introduce a generalization 
of the classical theory of sampling of d-dimensional fields on 
countable sets of points to a theory of sampling on countable 
sets of K-dimensional manifolds on M'' where k < d. Some 
examples of sampling schemes for sampling a field in M.^ are 
illustrated in Figure II-AI Figure |2(a)| depicts sampling on a 
lattice. Figure |2(b) sampling on a set of equispaced parallel 
lines, and Figure 2(c)| sampling on a set of equispaced parallel 
planes, corresponding to k = 0, k = 1, and k = 2 respectively. 

Initially we focus on the k = 1 case since our primary 
motivation comes from the problem of mobile sensing. The 
trajectory of a mobile sensor can be interpreted as a 1- 
dimensional manifold, or in other words, a curve, through 
space. Mobile sensing has an advantage over classical static 
sensing in that a single sensor can be used to take measure- 
ments at several positions within an area of interest ifTOl . 
Moreover, in some applications |11| moving sensors can 
sample the fields along their paths at high spatial frequencies 
thereby reducing the amount of spatial aliasing introduced in 
the samples. Furthermore, as we point out in lfT2l and lfT3l . 
a moving sensor admits filtering over space in the direction 
of motion of the sensor whereas no such spatial filtering is 
possible in the case of static sampling. Such spatial filtering 
helps in reducing the amount of aliasing and the contribution 
of out-of-band noise in the reconstructed field. 



B. Main results 

Our results for k = 1 provide guidance for designing 
trajectories for mobile sensors moving through space. Fo- 
cusing initially on straight line trajectories we identify some 
configurations of straight line trajectories that admit perfect 
reconstruction of bandlimited fields. This can be interpreted as 
a generalization of the classical Nyquist sampling criterion for 
sampling on a lattice to sampling on lines. We also formulate 
the problem of designing sampling trajectories with minimal 
path density that admit perfect reconstruction of bandlimited 
fields and obtain partial solutions to this problem, restricting 



ourselves to specific classes of straight line trajectories. Our 
main results for straight line trajectories are the following: 

• Necessary and sufficient conditions on a union of < 2 
sets of equispaced parallel lines that admit perfect recon- 
structions of bandlimited fields in (see Theorem 12.41 ). 

• Optimality of a single set of equispaced parallel lines 
from among unions of sets of equispaced parallel lines 
for sampling bandlimited fields in (see Theorem l3.2l ). 

• Optimum configuration of a uniform set of parallel lines 
for sampling bandlimited fields in R'^ (see Corollary 

EH). 

• Optimality of a uniform set of parallel lines from among 
all trajectory sets that visit all points in a sampling lattice 
(see Propositions 13.31 and I3.5l l. 

We also consider non-affine trajectories such as concentric 
circles and interleaved spirals and discuss some known results 
and some new results on sufficient conditions for perfect 
reconstruction. These results are based on an application of 
Beurling's theorem on sufficient conditions for non-uniform 
sampling. 

We then consider higher dimensional sampling manifolds 
corresponding to k > 1. We introduce a manifold density 
metric for sampling manifolds that generalizes the path density 
metric for sampling trajectories. We generalize our results on 
sampling on lines in M."^ to sampling on hyperplanes in R'' 
which are K-dimensional manifolds where k = d — 1. Our 
main results for hyperplanes are: 

• Necessary and sufficient conditions on a union of N < d 
sets of equispaced parallel hyperplanes that admit perfect 
reconstructions of bandlimited fields in M.'^ (see Theorem 

• Optimality of a single set of equispaced parallel hyper- 
planes from among unions of sets of equispaced parallel 
hyperplanes for sampling bandlimited fields in M.'^ (see 
Theorem |4~4]|. 



The dependencies between various results is summarized in 
the form of a table in Table U in the conclusion section of the 
paper 



C. Related work 

Although there has not been any past work specifically 
on the problem of designing optimal sampling manifolds or 
trajectories, some reconstruction schemes based on measure- 
ments taken along concentric circular trajectories have been 
proposed by Tewfik et al. 1141 and Myridis et al. ITSll . Various 
samphng trajectories have also been studied in the context 
of Magnetic Resonance Imaging (see, e.g., 121, lT6l . ITtI . 
ITSl ). However, the literature on trajectory design for MRI 
is of a different flavor from ours, since the primary focus 
in these works is to suppress noise in the reconstruction and 
they do not aim for exact reconstruction. To the best of our 
knowledge this paper is the first to introduce the notion of 
path density and manifold density and the notions of optimal 
sampling trajectories and optimal sampling manifolds for exact 
reconstruction of bandlimited fields. 
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D. Notations and conventions 

Most of the notations and conventions we use are described 
where they first arise. In Appendix|A]we provide a detailed list 
of symbols. Here we present some commonly used notations. 

We denote a field in d-dimensional space by a complex- 
valued mapping / : R'' i-^ C. For a field /(.), we define its 
Fourier transform F as 



f(r) exp(— r))dr, 



uj e 



(1) 



where i denotes the imaginary unit, and {u, v) denotes the 
scalar product between vectors u and v in Mf^. We use Bn to 
denote the collection of fields with finite energy such that the 
Fourier transform F of / is supported on a set O C K'^, i.e.. 



Bn {/ e L^R'^) : F{lu) = for <^ n}. 



(2) 



For f2 C K'' we use 51 to denote its interior in M''. Also for 
any s G M'' we use il{s) to denote the set obtained by shifting 
by s, defined as 

n{s) ■={xeR'^ -.x- se n}. (3) 

We use the following notations for order statistics. When 
o — >■ CX3 we say that a real-valued function h{a) G 0(a) if 
there exists fc, A e R such that for all a > A, we have, 

|/i(a)| < ka. 

Similarly, we say h{a) G o(a) if for every e > there exists 
A g R such that for all a > ^, we have, 

|/i(a)| < ea. 

We use S{.) to denote the Dirac-delta function in d- 
dimensions. The value of d is understood from context. We 
use 5ij to denote the Kronecker delta function. For vectors 
v,w £ M'', we use ||u|| to denote the Euclidean norm of v and 
{v,w) to denote the Euclidean inner product. 

E. Outline 

The rest of this paper is organized as follows. In Section 
|II]we formally define the notion of sampling trajectories and 
introduce the Nyquist criterion for sampling trajectories. We 
describe various designs of sampling trajectories and provide 
conditions on these trajectories for perfect reconstruction 
of bandlimited fields. We provide results on trajectory sets 
composed of straight lines, circles, and spirals. In Section Hill 
we study the problem of optimizing sampling trajectories in 
terms of the path density metric and present some optimal- 
ity results from certain restricted configurations of sampling 
trajectories. In Section |IV] we discuss sampling on higher- 
dimensional manifolds. We generalize our results on sampling 
on straight line trajectories to sampling on hyperplanes in 
high-dimensional spaces. We discuss reconstruction schemes 
in Section |V] and conclude in Section |VT] 

II. Sampling trajectories 

We now discuss the problem of designing trajectories for 
sampling. We consider natural generalizations of sampling 
lattices and unions of sampling lattices to trajectories. 



A. Preliminaries 

We use the following terminology in this section. A trajec- 
tory Pi in M'' refers to a curve in R''. We represent a trajectory 
by a continuous function p(.) of a real variable taking values 



P ■ 



A simple example of a trajectory is a straight line defined by 
p{t) = w + vt for some w,v G K'^. A trajectory set P is 
defined as a countable collection of trajectories: 



P = {p, : ^ e 1} 



(4) 



where I is a countable set of indices and for each i G 1, Pi 
is a trajectory in the trajectory set P. A simple example of a 
trajectory set is a countable collection of parallel lines through 

We introduce a natural generalization of the sampling 
density metric that characterizes the sampling efficiency of 
sampling lattices. In a sampling scheme using mobile sensors 
it is much more difficult to increase the spatial density of 
trajectories than to increase the sampling rate along the trajec- 
tories. Hence, unlike in classical sampling theory, the density 
of sampling points in space is no longer an appropriate metric 
for quantifying the efficiency of a mobile sampling scheme. A 
more reasonable metric is the total length of the trajectories 
required to span a field of given spatial volume. Let i?^ and 
B^{x) denote d-dimensional spherical balls of radii a centered 
respectively at the origin and at a point a; G R'^. For any given 
trajectory set P we denote its path density by £{P) defined 
as follows: 

sup^gRd VP{a,x) 



i{P) ;= lim sup ■ 



(5) 



Volrf(a) 

where V^{a,x) represents the total arc-length of trajectories 
from P located within the ball B^{x) and Void (a) represents 
the volume of the ball. Clearly, Vol2(a) = 7ra^ and Vol3(a) = 
|7ra^. For a trajectory set composed of differentiable functions 
Pi we note that T)^{a,x) can be explicitly calculated as 



/ 


dpi{t) 


lteTi(a,x) 


dt 



dt 



(6) 



where \\x\ 



represents the Euclidean norm of a; e 1 

T,{a,x) :={<eM: \\p,{t) - x\\ < a} 



and 



represents the portion of trajectory pi{.) that lies within B^{x). 

We say that a set of points A C R"^* is uniformly discrete 
if we have inf{||j: — y|| : x,y G A,x y} > 0, i.e., there 
exists r > such that for any two distinct points x,y G A 
we have ||x — > rQ For a set 51 C R'', let Bn denote 
the class of fields bandlimited to 51 as defined in Let 
Aq denote the collection of all uniformly discrete sets A C 
R^ which have the property that any field f E Bn can be 
reconstructed exactly from its values on A, i.e. f E Bn is 
uniquely determined from {f{x) : x £ A}. Classical sampling 
theory is primarily concerned with the elements of An, e.g., 

'For example lattices in R'' are uniformly discrete, but a sequence in R'' 
converging to a point in R"* is not. 
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Nyquist sampling lattices (T9\ . We now introduce the desirable 
properties of sampling trajectory sets. 

Definition 2.1: A trajectory set P of the form ^ is called 
a Nyquist trajectory set for 17 C M"* if it satisfies the following 
conditions: 

(CI) [Nyquist] There exists a uniformly discrete 
collection A of points on the trajectories in P 
such that A admits perfect reconstruction of 
fields in Bn, i.e., A C te(t) : i e I,t € M} 
and A G Aq. 

(C2) [Non-degeneracy] For any x E W^, there is 
a continuous curve of length no more than 
'D^{a,x) + o{a'^) that contains the portion 
of the trajectory set P that is located within 

We also introduce a special notation for the collection of all 
Nyquist trajectory sets: 

Definition 2.2: We use JVn to denote the collection of all 
Nyquist trajectory sets for fl, i.e., Afn is the collection of all 
trajectory sets P of the form (|4| that satisfy conditions (CjT) 
and (C©. 

The condition (Cjl) ensures that the entire field can be 
reconstructed exactly from samples taken on the trajectories. 
This is the direct analogue of the Nyquist condition for 
sampling on points. The restriction that only uniformly discrete 
collections of samples are allowed is a standard one in non- 
uniform sampling theory (see, e.g., Il20ll . ifTSl ). There is a 
subtle reason for adopting this restriction in this work. We 
know that a bandlimited signal is an entire function and hence 
can be completely recovered (see, e.g., 11211 . Il22l ) using its 
values on an interval or using its values on a convergent 
sequence of points. Such results can be used to construct a 
trajectory in K.'' with finite total length such that bandlimited 
fields can be perfectly reconstructed using the field values on 
the trajectory. However, in reality it is impossible to accurately 
measure the field values on an interval or on a convergent 
sequence of points. The condition (C|2) ensures that the path 
density metric does indeed capture the total length that needs 
to be traversed by a single moving sensor using the trajectories 
in P for sampling. This avoids degenerate static cases like the 
situation in which every trajectory pi corresponds to a single 
point in a sampling lattice for the field /. Such a degenerate 
trajectory set has path density equal to zero, but it is not 
possible to visit all of these points by using a single sensor 
In fact, we know from Nyquist sampling theory that if we 
had 0{a'^) sensors available for sampling spherical regions 
of radius a, it may be possible to sample bandlimited fields 
without any movement at all. 

In a practical deployment, it is not possible to take mea- 
surements of the field at all points along a continuous path 
because a continuous path has an infinite number of points. 
However, if the sensor moving along a trajectory is exposed to 
a bandlimited function of time it is possible to reconstruct the 
entire field along its path from uniformly spaced samples. This 
motivates the following additional desirable condition (K[l) of 
a trajectory set P : 




Fig. 3. Illustration of the result of Lemma ITT] for fields in R^. The vector 
V represents the velocity vector in the spatial domain and p and p represents 
the limits of the spectrum of /(p(.)) in the frequency domain. 




(a) Uniform set (b) Union of uniform (c) Concentric circles 
sets 



Fig. 4. Three choices of sampling trajectory sets for : uniform sets, union 
of uniform sets, and concentric equispaced circular trajectories. 

(Kl) For any field / e Bn the 1-dimensional signal 
f{pi{-)) is bandlimitecQ for all pi G P. 
The following lemma shows that trajectory sets composed of 
straight lines satisfy condition (K[T]|. 

Lemma 2.1: For w,v E M'* let p{t) = w + vt denote a 
straight line trajectory parameterized by t. Then for / e Bn 
the function f{p{.)) is bandlimited to the set [p,J>] where 

p = min{(a;,w) : a; € 12} and p ~ ma,x{{uj,v) : uj G f2}. 

□ 

We do not prove this lemma since it follows easily from 
elementary properties of the Fourier transform (see, e.g., |[23ll ). 
The hmits of the spectral support p and p are illustrated in 
Figure [3] Because of this desirable bandlimitedness property 
of f{p{-)) for straight line trajectories, almost all the trajectory 
sets that we study in this paper are collections of straight line 
trajectories. 

B. Sampling trajectories for 

We now present some simple examples of trajectory sets on 
and compute their path densities. We discuss conditions 
required to ensure that these trajectory sets are Nyquist trajec- 
tory sets for specific choices of 17 C M^. These results can be 
interpreted as a generalization of known results on conditions 
on sampling lattices lfT9ll and unions of sampling lattices ll24ll 
for perfect reconstruction of bandlimited fields. 

1 ) Uniform set: The most natural choice of a trajectory set 
in M? is a collection of equispaced straight line trajectories of 
the form 

P^{p,:ie1} with pi{t) = w + iAv-'^ vt (7) 

-i.e. has compact support in the Fourier domain 
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where A G R+, and w,v are fixed vectors in witii 
denoting the unit vector orthogonal to v. These trajectories 
are lines oriented parallel to v and are spaced A units apart. 
We refer to such a set of equispaced parallel line trajectories as 
uniform sets in R^. Uniform sets form a periodic configuration 
of straight Unes and hence can be interpreted as a natural 
generalization of sampling lattices to sampling trajectory sets. 
It is also immediate from Lemma IZTI that uniform sets satisfy 
the desirable condition (K[l). 

As an example let P be a uniform set parallel to the x-axis 
with the spacing between adjacent lines equal to A, defined 
by P = {pi : i E Z} where 

p^{t) = {t,Aif,ieZ,teR. (8) 



Such a uniform set is shown in Figure |4(a)| Suppose is a 
compact convex subset of M^. Let 



A* := sup{A ; fin O 




is empty} 



where we use the convention of (|3]l for a shifted version of a 
set. It can be easily shown via classical sampling theory |fT9l 
that P forms a Nyquist trajectory set for O if A < A*. More 
precisely, it can be shown that 



and 



P e A/'o if A < A* 



P ^ A/h if A > A* 



(9) 



(10) 



These results, as well as analogous results for general uniform 
sets, follow as special cases of Theorem 12.41 which we prove 
later in the paper. We now illustrate the main idea behind the 
result via a simple visual proof. Since the sensors can measure 
the field at high resolutions, we can assume that for any e > 
the sensors can measure the field value at points of the form 
Pi{je),j e Z. In this case, we have access to field samples 
f{{je,iA)^),i,j £ Z, which correspond to samples on a 
rectangular lattice. We know from classical sampling theory 
that the sampled impulse stream. 



has a Fourier spectrum Fs composed of spectral repetitions of 
F over a reciprocal lattice of points in the Fourier domain: 

V- / / 27r 27r^ 
P,H = ^P c.- .7-,*- 



Such a spectrum with repetitions is illustrated in Figure |5(a)| 
where fl is chosen to be a circle. Now, e can be made 
arbitrarily small. In the limit as e — >■ 0, the sampled spectrum 
satisfies 




This means that the spectral repetitions along the x-axis no 
longer overlap, leading to Figure [5(b)| Now perfect recovery is 
possible whenever the spectral repetitions along the y-direction 
do not overlap, or equivalently when A < A*. 




(a) Sampled spectrum 



(b) Limit as e 



■ 



Fig. 5. Sampled spectrum of a two-dimensional field bandlimited to a circle 
S7. The support of the sampled spectrum is composed of periodic repetitions 
of n on the plane. 



It is straightforward to generalize the results of (|9]l and dTol ) 
to general uniform sets of the form (|7]i by interpreting these 
sets as shifted and rotated versions of the simple example 
we considered in The exact conditions can be explicitly 
obtained from Theorem 12 . 41 which is stated in the next section. 
In the following lemma we characterize the path density of a 
uniform set. 

Lemma 2.2: The path density of a uniform set P of the 
form in (|7]i is given by 



Furthermore, P satisfies condition (C©. 



(12) 



□ 



We provide a proof in Appendix |B] 

The analyses for uniform sets can be extended to a more 
general configuration of straight line trajectories composed of 
unions of uniform sets which we discuss next. 



2) Union of uniform sets: For vectors Wi 



Vi e 



let P 



denote the uniform set defined by P; = {pi.j : j Eli} where 

p,,j(t) = jA^vl ^tv^,j e Z,f e M, (13) 

where is a unit vector in M? orthogonal to v and A^ > 0. 
Thus Pi is a uniform set oriented parallel to Vi with the spacing 
between adjacent lines equal to A^ apart. The vector Wi is just 
an offset from the origin. We define the vector Ui E as 

u. = (14) 

If we assume that we are given samples of the field / at points 
of the form pi,j{ke), j, k € Z, then the sampled spectrum 
from these samples corresponds to repetitions of the field 
spectrum P. As in the example we considered in Section lTl-Bll 
if we choose a small enough value of e, then it can be shown 
via classical sampling results |fT9l that the sampled spectrum 
comprises spectral repetitions in one direction. We obtain the 
following spectrum analogous to ( fTTl i: 



iez 



(15) 
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Now suppose we have N such uniform sets {Pi 1 < i < 
N} of the form ( fT3] ). Let P denote the union of all the sets: 

N 

P:^[jn- (16) 

i=l 

An example of such a trajectory set for iV = 2 is depicted in 
Figure |4(b)| It follows from Lemma 12.21 that the path density 
of this trajectory set is given by 

i—l * 

We now seek the conditions under which P forms a Nyquist 
trajectory set for il. Since we now have samples from N 
different uniform sets, it may be possible to reconstruct the 
field although the sampled spectrum from each individual 
uniform set given by ( fTSl l is aliased. Let Q C denote the 
set of points 

Q ■■= |E(-i)''y ■■ ^» e 1}, 1 < » < • (18) 

The following proposition gives a necessary condition that 
must be satisfied by any Nyquist trajectory set for Vl that is a 
union of uniform sets represented in the form ( fT6b . 

Proposition 2.3: Let il C be a compact convex set. 

o 

Suppose Q C r2(.s) for some s e R^. Then P ^ Mn. □ 

o o 

Here O denotes the interior of O and Q,{s) denotes the shifted 
version of the interior as defined in (|3]l. The above result can 
be proved by constructing a sinusoidal field that vanishes on 
all the lines in P, and has a Fourier transform supported on 

o 

C Vt. We provide a proof in Appendix ICl The result 
is equivalent to the fact that every P e Mn must satisfy 

Q ^ h.{s) for all s e . 

This necessary condition must be satisfied by all unions of 
uniform sets P that form Nyquist trajectory sets for 51. When 
iV < 2, it can be shown that the tightest necessary condition 
given by Proposition 12.31 is also sufficient. Our main result 
in two dimensions is the following theorem which provides 
necessary and sufficient conditions to ensure that P forms a 
Nyquist trajectory set for SI when N < 2. 

Theorem 2.4: Let C be a compact convex set. Let 
P denote a union of uniform sets expressed in the form of 
([T&t with = 1 or 2. For N — 2 assume that vi and V2 are 
non-collinear. Then we have 

P eAfnif ^(s) for all s e (19) 

and, 

P ^ Ml if Q C n{s) for some s e R^ (20) 

where Q is defined in (fTsl l. □ 

The result of ( |20l i clearly follows from Proposition 12.31 The 
rest of the theorem is proved in Appendix |D] 

Theorem 12.41 gives us necessary and sufficient conditions 
on unions of uniform sets that admit perfect reconstruction. 



The conclusion shows that it is possible to recover the field 
completely using readings on multiple uniform sets even when 
the measurements on each individual uniform set may be 
aliased. The result of Theorem l2.4l can be generalized to higher 
dimensions where sampling on lines is replaced with sampling 
on hyperplanes. This generalization is contained in Theorem 

The results of Theorem 12.41 can be interpreted as special 
cases of classical sampling on a discrete collection of points. 
Since uniform sets are composed of straight line trajectories 
we know from Lemma IZTI that the restriction of 
the bandlimited field / to any line pij in the uniform set 
of ( fT3] l is bandlimited. Clearly, for a fixed i, the bandwidths 
of the f{pi.j{.)) are identical for all j. Now suppose that 
the sensors moving along each of the lines in the set Pi 
take uniform spatial samples that are apart. It follows 
via the bandlimitedness of f{pi,j{.)) that the field values 
{f{pi,j{t)) : t £ M.} can be recovered exactly from the 
samples {f{pij{mei)) : m,j E Z} provided e,; is small 
enough. It is also clear that the points {pi.j{mei) : m,j £ Z} 
lie on a shifted version of a rectangular lattice. Thus for 

= 1, the conditions given in Theorem 123] can be interpreted 
as the condition for perfect recovery under a special case 
of sampling on a shifted rectangular lattice, viz., when the 
sampling interval along one direction of the lattice is arbitrarily 
small. We have already seen this interpretation in Section 
HLBT] where we illustrated the idea using Figure III-Bll We 
can extend this interpretation to the N ^ 2 case. For N = 2, 
the collection of sample locations from the two uniform sets, 
viz., the collection {pij{mei) : m, j G Z, 1 < i < 2} defines 
a union of two shifted rectangular lattices. Hence the result 
of Theorem 12.41 can be interpreted ^ as the condition for 
perfect recovery under a special case of sampling on a union 
of two shifted lattices 1241 . when the sampling intervals are 
arbitrarily small. We note that in this example, the samples 
taken on each individual shifted lattice are aliased but given 
all sets of samples, perfect recovery is possible. The value of 
ei required can be determined from the bandwidth of f{pi,j{.)) 
via Lemma 12.11 

We also note that the result of Theorem 12.41 on unions of 
uniform sets is similar in spirit to known works on the related 
problem of sampling on the union of shifted versions of a 
lattice (see, e.g., (ES), ||26|, lEll). These works (e.g., E?] 
Example 2]) provide necessary and sufficient conditions for 
sampling on unions of shifted versions of a lattice, whereas, 
as we argued in the previous paragraph, our results provide 
conditions for sampling on certain unions of shifted non- 
identical lattices. Furthermore, as can be seen from the proof 
of Theorem 12.41 the condition 

Q ^ n{s), for all s S R^ 

is equivalent to the condition for invertibility of a linear system 
of equations relating the values of the Fourier transform of 
the field. This is analogous to the invertibility condition of 
Papoulis [27 1 . 

We conclude this section on unions of uniform sets with an 
example that illustrates the result of Theorem 12.41 for mutually 
orthogonal uniform sets under two different choices of ft. 



7 



V 



A 



(a) Orthogonal sets of trajectories. 





(b) Sampled spectra from the two sets. 



(c) Critical sampling. 



Fig. 6. Two sets of mutually orthogonal uniform sets and the sampled spectra from samples taken on these trajectories discussed in the first part of Example 



Example 2.1 (Mutually orthogonal uniform sets): 
Consider a trajectory set P composed of the union of 
two uniform sets, one oriented parallel to the x-axis spaced 
Ai units apart and another parallel to the y-axis spaced A2 
units apart. 

(i) Suppose Ai = A2 = A. Such a set of trajectories 
is illustrated in Figure [6(a)] with the trajectories parallel 
to X-axis (y-axis) colored in black (red). Suppose il 
is a circular disc of radius p centered at the origin. 
The sampled spectra corresponding to the samples from 
each set of trajectories is illustrated in Figure |6(b)| 
From Theorem 12.41 it follows that in order to guarantee 
perfect recovery we must have at least one point of 
the form {zt^,zt^)^ must lie outside H. Hence the 
maximum value of A required to ensure that the field 



can be perfectly reconstructed is A* 



V27r 



From 



simple geometry it can be seen that when A = A*, 
the first spectral repetitions along the x-axis touch the 
first spectral repetitions along the y-axis as illustrated in 
Figure [6(c)| Thus when A < A*, these first repetitions 
along the two axes do not overlap and thus every point 
uj £fl satisfies either F{uj) = F}{uj) or F{uj) = F^{uj) 
and thus the whole spectrum e ft can be 

recovered using both sampled spectra, 
(ii) If is not a circular disc, the conditions are a bit 
more complex to state. Suppose Ai = 2A2 = 2A. Such 
a set of trajectories is illustrated in Figure [7(a)| with the 
trajectories parallel to x-axis (y-axis) colored in black 
(red). Consider this example from f2, Example 3.2] when 
r2 has the shape of a right- triangle given by 51 := {a; G 
: Ljy > 0,\ujx\ + < p}- The sampled spectra 
corresponding to the samples from each set of trajectories 
is illustrated in Figure |7(b)| In this case Theorem 12.41 
predicts that perfect reconstruction is possible whenever 
A < A* := When A = A* the sampled spectra are 
as shown in Figure [7(c)| We have partitioned fl into seven 
distinct portions, each portion being characterized by the 
nature of spectral overlap in the two sampled spectra. We 
see that the portions of O such as and O7 in Figure 
|7(c)| are aliased in the spectra from both the uniform sets. 
Nevertheless, Theorem 12.41 shows that whenever A < 



A* := ^ perfect reconstruction is possible. In this case, 
it can be verified that the sampled spectra at the different 
portions of ^l satisfy the following relations 

F^iuj) = F{uj)+F{uj + ui), weflgUflr 
F^{uj) = F{uj), w e 111 U U fla 

where F^ and F^ denotes the sampled spectra from the 
trajectories parallel to the x-axis and y-axis respectively. 
Thus the original spectrum can be recovered as 

= (w)xoiUi22U03(w) +-Fi(t^)xn4UJ25(^) 

where xa{-) denotes the indicator function of set A. 
Hence the original field can be recovered by inverting 
the Fourier spectrum. □ 



The two cases considered in Example 12. II are distinctly dif- 
ferent. In the first case, we saw that the conditions for perfect 
recovery is the same as the requirement that every part of 
n is unaliased in at least one of the two spectra. Whereas, 
in the second case, we saw that perfect recovery is possible 
even when some portions of the spectrum are aliased in both 
sampled spectra. 

3) Concentric equispaced circular trajectories: We now 
consider an example of a set of non-affine trajectories. Suppose 
fi is a circular disc of radius p centered at the origin. Let 
C = {ci : i E NU {0}} where q denotes a circular trajectory 
of radius iA centered at the origin. Such a trajectory set is 
shown in Figure |4(c)| It is easily verified that condition (Cl2) 
is satisfied by this trajectory set. It is also known from ifTSll that 
any field bandlimited to ft is reconstructible from its values 
on these trajectories whenever A < However, this does 
not verify whether condition (CH) is satisfied because exact 
reconstruction of the field using the scheme in ifTSll requires 
the field values at all points on the continuous circular curve. 
Nevertheless, one can use the Beurling frame theorem to show 
that indeed condition (CHJl is also satisfied when A < ^. 
Below we give a version of Beurling's theorem for isotropic 
bandlimited fields in R^, i.e., fields bandlimited to a circular 
disc in the Fourier transform. This version of the theorem is 
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(a) Orthogonal sets of trajectories. 



(b) Sampled spectra from the two sets. 




Y 

(c) Critical sampling. 

Fig. 7. Two sets of mutually orthogonal uniform sets and the sampled spectra from samples taken on these trajectories discussed in the second part of 
Example 12.11 



taken from ifTSl . The slight difference from the statement in 
IfTSl is due to the different definition of the Fourier transform. 

Theorem 2.5 (Beurling's Covering Theorem): Let i7 C 
be a circular disc of radius p centered at the origin. For any 
uniformly discrete set A C M'* let 



^^(A) 



sup inf II x — y\\ 



If q{A) < then A e An , i.e., any / e i5n is uniquely 
determined from {f{x) : x E A}. □ 

We can now use this theorem to verify that condition (Cjl) is 
satisfied by C. For e > small enough we note that on each 
circle Ci we can choose a discrete collection of equispaced 
points such that the arc-length between nearest neighbors on 
the circle lies in the interval [|e, e]. Let A denote the set of 
all such points together with the origin. It follows that A is 
uniformly discrete and satisfies g{A) < i(A + e). Now if A < 
we can choose e small enough to ensure that g{A) < f^. 
Thus by Beurling's covering theorem it follows that C G Afn 
whenever A < 

We now compute the path density for this trajectory set. 
Clearly, the maximum density of trajectories is at the origin 
and hence the supremum in Q is achieved at the origin. A 



disc of radius a with center at the origin contains a total of 
A' « ^ concentric circles separated by a radial distance of A. 
This leads to a total trajectory length of TrAK{K + l) whence 
we get £{C) = ^. Since we require A < ^ to ensure that 
C G Afn it follows that the lowest possible path density among 
trajectory sets that satisfy this condition is achieved when A 
meets the upper bound and is given by ^. 

We also note that this trajectory set does not satisfy condi- 
tion (KU) in Section Hi- Al because the signals along the circles 
are not bandlimited. Nevertheless, these signals have a finite 
essential bandwidth as established in Lemma IZ6] below. 

Lemma 2.6: Let / G Bn where £7 C M^. Consider a sensor 
moving at a constant angular velocity of v along a circle 
of radius a centered at the origin. The time-domain signal 
= f {{a cos{vt) , a sni{vt))) is essentially bandlimited to 
[-ps^Ps] where ps = v{l + asup^g^^ \\uj\\). □ 

We provide a proof in Appendix |E] 

We note that the signal s{t) in Lemma l276l is periodic in 
time. Hence it can be represented by a Fourier series expansion 
s{t) ~ YlV=-oo cxp(iA:j/<). The coefficients Sk can be 
computed explicitly using Bessel functions as: 

Sk = - [ F{uj)exp{\kl3)Jk{a)doj (21) 
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where a ~ o.\\^\\ and (3 ~ tan 
coefficients are negUgible for \k\ > k 
s{t) can be approximated as 



1 'pjjg pourier series 
\^~\ . Thus the signal 



s{t) 



k 

E 

k=-k 



Sk exjp{\ki't). 



(22) 



Hence if we have measurements of s{t) at a finitely many 
points of time, then we can use the approximation of ( l22b to 
estimate the periodic s{t) from these samples. 

4) Union of spiral trajectories: Another type of non-affine 
trajectories are spirals. The geometry of spiral trajectories are 
particularly convenient to use in MRI applications and hence 
they have been studied extensively in the field of MRI (see, 
e.g, 1281 . ifTSll ). For example the work ifTSl contains some 
results on perfect reconstruction of bandlimited fields from 
their measurements taken on spiral trajectories and unions on 
spiral trajectories. Suppose ft is the circular disc of radius p 
centered at the origin. Consider a trajectory set composed of 

interleaved spirals of the form 



5P = {sK :^G {0,l,...,iV-l}} 
where spi are Archimedean spirals of the form 

spt{t) 



(23) 



cicos(27r(< - i/N)) 
cts\\i{2TT{t- i/N)) 



for i e {0, 1, . . . , — 1}. An example of such a trajectory set 
for = 3 is shown in Figure |8] It follows via ifTSl Example 
3] that SP satisfies condition (Cjl) whenever < ^- This 
result is established by applying Beurling's covering theorem 
(Theorem l2.5b to identify a uniformly discrete set A of points 
on the spirals satisfying A g An- Moreover, since SP is 
composed of a finite number of continuous trajectories it 
follows trivially that the condition (C|2) is satisfied by this 
trajectory set. Hence it follows that SP E Afn whenever 

C ^ TT 

N ^ p- 

In the following lemma we compute the path density of the 
trajectory set of interleaved spirals. 

Lemma 2.7: The path density of the union of spiral trajec- 
tories of (O satisfies e{SP) = f . □ 

We provide a proof in Appendix |F] Since we require jj' < ^ 
to ensure that SP E Afn it follows that the lowest possible 
path density among trajectory sets that satisfy this condition 
is achieved when c meets the upper bound and is given by ^ 
for all values of N. 

Before we conclude this section, we mention a caveat on 
the use of Beurling's covering theorem (Theorem 12.5b for 
identifying optimal spacings in trajectory sets. This theorem 
gives only a sufficient condition on sampling sets. Although 
we used the conditions for this theorem to find conditions 
on the trajectory sets of Sections III-B3I and III-B4I required 
to ensure that they are Nyquist trajectory sets, these are only 
sufficient conditions. Some trajectory sets that do not satisfy 
these conditions may also be Nyquist trajectory sets. For 
instance, in Example 12.1 1 (i) of Section iri-B2l we saw that the 
maximum possible spacing possible for a union of equispaced 
orthogonal uniform sets is A* = However, for a union of 




Fig. 8. A trajectory set in 
Arcliiniedean spirals. 



composed of a union of three interleaved 



equispaced orthogonal uniform sets satisfying the conditions 
of Beurling's theorem we need to have A < ^. This proves 
the existence of Nyquist trajectory sets that do not satisfy the 
conditions of Beurling's theorem. 



C. Sampling trajectories for where d > 3 

For (i > 3, we consider only trajectory sets composed of 
periodically spaced parallel straight lines analogous to the 
uniform sets we considered in Section III-Bll for d = 2. 
Let {vi,V2, ■ ■ ■ ,Vd} denote a basis for M."^ such that Vd 
is a unit vector orthogonal to the hyperplane spanned by 
{vi,V2, . . . , Wd-i}- In other words {vi,Vd} = S.i^ for all i < d. 
Consider trajectories pm of the form 



i=l 



tVd, 



t e 



(24) 



where m = (mi, m2, . . . , md-i)'^ E Z'' ^. Let P denote the 
trajectory set 

P = {p„ :meZ'^-i}. (25) 

We refer to such a set of periodically spaced parallel line 
trajectories as uniform iefjj in R''. A simple example of a 
uniform set in M'' is illustrated in Figure |9] corresponding to 



Vl = 







( 


:) 


,V2 = 


A2 






I 



1'3 = 




(26) 



We use Fn C Nn to denote the collection of all uniform 
sets that form Nyquist trajectory sets for Q,. In other words, 
Fn is the collection of uniform sets in R'^ that satisfy 

'For an exact generalization of uniform sets in M^, we have to include an 
additive shift of some vector w G to Pm{t)- We avoid it here to keep the 
presentation simple. Generalization of the results to non-zero shifts is trivial. 
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Fig. 9. A simple example of a uniform set in ' 



conditions (C[T]| and (C©. The following theorem provides 
conditions on the vectors {wi, U2, . . . , Wd} and on the set $7 so 
that P e TVh- 

Theorem 2.8: Let P denote the uniform set defined in dZST l. 
Let {ui, . . . , Md-i} denote vectors in R'' satisfying {ui, Vj) = 
27r(5,y for 1 < I < d - 1 and 1 < j < d and let C R'^ 
denotes a compact convex set with a point of symmetry at the 
origin. Then we have P G Afu if 



i ^iUi ^ ri, for all m e 



7(i-l 



\{o}. 



(27) 



The condition dZTj ) is also necessary in the sense that if there 

exists m e Z"^'^^ \ {0} such that i "^i'^i £ where 17 

denotes the interior of fl in M'', then P ^ A/j^. □ 

We provide a proof in Appendix |G] The crux of the proof lies 
in the fact that the sufficient condition to ensure that P satisfies 
condition (Cjl) is that for e small enough the lattice of points 
defined by {X^it^i m-iVi + mdevd : m G Z"^} forms a sampling 
lattice for ft, i.e., admits perfect reconstruction |fT9l of fields 
in Bn from its samples on the lattice. Since ft is convex and 
symmetric about the origin this condition is equivalent to dZTl i. 
We illustrate the result for a simple example below. 

Example 2.2 (Rectangular uniform sets): Consider the uni- 
form set shown in Figure |9] corresponding to the basis vectors 
of (|26] |. In this case it is easily verified that 



Ml 



27r 
Ai 






, U2 




(28) 



Hence, for this trajectory set, the condition in ( |27] | is equivalent 
to the condition that the vectors 



Ai 




and 





A2 



generate a sampling lattice for the set Vl^y, the intercept of VI 
with the xy plane defined as VL^y = {(•'J;,?/) : (a;, y,0) € fi}. 
It is easily verified from classical sampling results |fT9l that 
for compact convex sets VL with a point of symmetry this is 
exactly the condition required to ensure that there is no aliasing 
in the samples of the field measured at points of the form 
{mivi + m2V2 + ETOaWs ; m £ 1?} for e small enough. Now 



suppose that $7 is a spherical ball of radius p centered at the 
origin. Then VL^y is just a circular disc so that the condition 
to ensure perfect reconstruction becomes 



max{Ai, A2} < -. 

P 



□ 



The path density of uniform sets in is characterized in 
the following lemma. 

Lemma 2.9: The path density of the trajectory set P defined 
in dZSl l is given by 

i{P) = ^- 

|det(G)|i 

where det(G) denotes the determinant of the c? — 1 x d — 1 
Gram matrix defined by Gij ~ {vi,Vj) ,1 < i, j < d — 1. □ 

We provide a proof outline in Appendix |H] 

As we did with Theorem l2.4l the results of Theorem l2.8l can 
also be interpreted as a special case of sampling on lattices. 
Since uniform sets in are composed of straight line trajec- 
tories we know from Lemma [ZTI that the restriction f{pm{ )) 
of the bandlimited field / to any line in the uniform set is 
bandlimited. It follows that the field along each line can be 
perfectly reconstructed using uniform sampling along the line, 
provided the sampling interval e is small enough. If we further 
assume that the sampling locations on the lines are aligned 
with each other, then the collection of the sampling locations 
on all lines can be expressed as {pm{nt) : n G Z, m e Z'^^^} 
where e is small enough. Clearly, this collection of points 
forms a sampling lattice in W^. Thus the results of Theorem l2.8l 
can be interpreted as the conditions for perfect reconstruction 
under a special case of sampling on a lattice in when the 
sampling interval along one direction is sufficiently small. The 
value of e required can be determined from the bandwidth of 
f{p„i{.)) via Lemma ITT] 

III. Optimal sampling trajectories 

As we argued earlier, the path density of a trajectory set 
captures the total distance required to be traveled per unit area 
for sampling spatial fields using a mobile sensor moving on 
the trajectory set. Hence it is of interest to characterize the 
optimal trajectory set for sampling fields that are bandlimited 
to a given set f2 C M''. We seek a solution to the following 
problem: 

min i(P). (29) 

In this section we identify partial solutions to the problem, 
solving it exactly for trajectory sets restricted to some subsets 
of AAn. 

We introduce some notation. Suppose {61, 62, . . . , fed} 
forms a basis for M''. A lattice b generated by the vectors 
{fei, 62, ■ • ■ , bd} is a collection of points in K'' of the form 




(30) 



For any lattice 6, it is known that (see, e.g., ||29l p. 276]) it is 
always possible to define a basis {ci, C2, . . . , Cd} of R'^ such 
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that b is generated by {ci, C2, . . . , Cd} and the basis vector ci 
is a vector with shortest length in b, i.e., 

ci € argmin ||a;||. 
xeb\{o} 

A lattice b is called a sampling lattice for a set 51 C if 
every bandlimited field g : M'' n- R with Fourier transform 
supported on 17 can be recovered perfectly using only the 
values of the field g at points in the lattice. The sampling 
density of a sampling lattice of the form (|30] | is defined as the 
average number of points per unit volume in M''. 

In the following lemma we present a simple result that 
identifies the shortest strategy for visiting all points in a 
sampling lattice for f2 C K''. 

Lemma 3.1: Let 6 be a lattice of points in M'' of the form 



(|30] | such that hi £ argmin^gj^|-Q-j. is a shortest vector of 
the lattice. Let Af'^ denote the set of trajectory sets P of the 
form dUi that visits all points in b: 



:= {P -.bd {p,{t) : i G I, t G M} and 
P satisfies condition (Cl2)}. 

Then the minimum in the following problem 

min e(P) 



(31) 



is achieved by the uniform set given by P = {pm '■ fn € 
Z''-!} where 

Pm{t) = ^ niibi+i + tbi,t e M. 
1=1 

Proof: Let P e and x £ K'^ be arbitrary. Since P 
satisfies condition (C|2) it follows that there is a continuous 
curve q of length I?^(a, x) +o(a'') that contains all the lattice 
points in 6 n B^{x). If one follows the curve q starting at 
one of its end-points, one will eventually visit all points in 
bn B^{x) before reaching the end of q. By construction the 
shortest path connecting any two points in & n B^{x) has a 
length no less than Hence the total length of q should 

satisfy 

Length(g)> (#(6ni?,^(x))-l)||6i|| 
where # denotes cardinality. This further implies that 

(#(6nBf(x))-l)||6i|j 



£{P) > lim sup 



Void (a) 



By the definition of P it is easily verified that choosing P = 
P achieves equality in the above relation. The desired result 
follows. ■ 
An immediate corollary of this lemma is the fact that for 
any set a uniform set has the shortest path density among 
all trajectory sets containing a sampling lattice for fl. In this 
section we use this lemma together with results from Section 
|II]to establish some optimality properties of uniform sets. 



A. Optimality for 

It is difficult to solve ( l29b exactly because it is difficult to 
characterize all the trajectory sets that satisfy conditions (CUll 
and (Cl2). However, as we show below, it is possible to identify 
the optimal trajectory set among those that can be written as 
a finite union of uniform sets like in (fTSI l. Such trajectory 
sets have the added advantage that they satisfy the desirable 
property of (K[i1i, as proved in Lemma IZTI 

Let £q C Afn denote the collection of trajectory sets Q in 
A/o such that Q is a finite union of uniform sets of the form 



Q^[jQi 



where for each i, Qi is a uniform set. We need the following 
definitions. For a nonempty compact convex set C M'' and 
any u let denote the distance between the two 

parallel supporting hyperplanes of f2 that are perpendicular to 
the vector u. We refer to B^(Jl) as the breadth of Vl in the 
direction u. The width of is defined by the relation 



min 



(32) 



A chord of is defined as the nonempty intersection of with 
a Une in M!^. For u G M'', yV"(ri) is defined as the maximum 
length of a chord of Vl parallel to u. The width >V(fi) can 
alternately also be interpreted as (see, e.g., Il30l ) 



W{n) = min Win). 



(33) 



In the following theorem, we identify the union of uniform 
sets with minimal path density such that the trajectory set is 
a Nyquist trajectory set for Q.. 

Theorem 3.2: For any nonempty compact convex set Vl C 
R^, let u be the minimizer in ( |32] |. and for e > let denote 
a uniform set given by = {p^ : j S Z} where 

Y 27r \ u 
where G M.^ is orthogonal to ii. Then G £n for all 



u-^tA G 



j G Z 



ee (0, 



2rr 



y^fjjYj) ™d is optimal in path density as e 

W(0) 



0, i.e.. 



lim eiP" 

e-)-0 



inf e{Q) 



2n 



□ 



The width >V(r2) and the optimizer u appearing in the state- 
ment of the theorem are illustrated in Figure pT)(a)| This result 
is established by identifying the union of uniform sets with 
minimal path density that satisfies the conditions of Theorem 
12.41 We provide a proof for the theorem in Appendix U This 
optimality result can be generalized to the case of sampling 
d-dimensional fields on unions of uniform sets of hyperplanes. 
This result is the topic of Theorem 14.41 

In short, Theorem l3.2l establishes the optimality of a uniform 
set P from all trajectory sets in E^. In particular, we have 
identified a sequence of uniform sets of trajectories indexed 
by e with path densities converging to the infimum path 
density from £q as e tends to zero. The limiting uniform set 
corresponding to critical sampling is obtained by plugging in 
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(a) Width VV(Sl) of a convex set f2 (b) Optimal uniform set for f2 corre- 
and the direction it along which Q sponding to critical sampling obtained 
is narrowest. The orientation u is the by plugging in e = in the solution 
minimizer in 132) . of Theorem 13.21 

Fig. 10. The width W{Q) of a convex set Q, and the limiting optimal 
uniform set corresponding to critical sampling given by Theorem |3.2| The 
optimal uniform set is oriented orthogonal to the direction u in which f2 is 
nan'owest. 



the best known path densities from the examples in Sections 
III-B3I and III-B4I also match the best path density of a uniform 
set when is a circular disc. However, it is not easy to 
verify whether this is indeed the optimal value of ( [29] l. We 
recall that in the classical case of sampling fields in Bn on a 
uniformly discrete collection of points, the best known lower 
bound on the sampling density is obtained from Landau's 
necessary conditions ISTI (see also 1201 Corollary 1]). It is also 
known that one can identify a collection of sampling points 
with a sampling density that is arbitrarily close to this bound 
13211 . l33l Cor 4.5]. It may be possible to extend Landau's 
conditions and the constructions of 1321 and l33l Cor 4.5] 
to Nyquist trajectory sets and thus obtain a solution to the 
problem ( l29l l of identifying the Nyquist trajectory set with 
minimal path density. 

In the following section we generalize some of the above 
optimality results for sampling trajectories in two dimensions 
to higher dimensions. 



e = 0. This critically sampled uniform set is a set of parallel 
lines oriented parallel to ii and spaced apart, as shown 

in Figure [TO(b)l 

In addition to the optimality property of Theorem 13.21 uni- 
form sets satisfy a different optimality property. In particular, 
uniform sets are optimal among all trajectory sets in Afn that 
contain all the points on a sampling lattice for ft, as the 
following result shows. 

Proposition 3.3: For any set 51 C K^, let Q e A/'o be a 
trajectory set that visits all points in a sampling lattice for f2. 
Then there exists a uniform set P e Mq such that £{P) < 

m- 

Proof: Let b denote the sampling lattice for Vl that is vis- 
ited by Q. It is immediate from Lemma ITTl that the minimum 
path density in the collection of all trajectory sets that visit h, 
is achieved by a uniform set. This uniform set automatically 
satisfies condition (CjT) since it contains a sampling set for Vl. 
Since uniform sets further satisfy condition (C© by lemma 
12.21 the result follows. ■ 
We now revisit the original problem ( l29l ) that we wished to 
solve. The results of Theorem l3.2l and Proposition 13 . 31 suggest 
that the uniform set given in Theorem 13.21 is optimal from a 
wide class of trajectory sets. It is tempting to consider the 
possibility that this optimality extends to a wider class so 
that this uniform set also solves ( l29l l. Some further evidence 
for this is obtained by identifying the trajectory sets with 
minimal path densities among those studied in Sections III-B3I 
and III-B41 Assume that Q, a circular disc Q, of radius p. In this 
case, we know from Theorem 13.21 that the best path density 
possible with a uniform set is ^. We also see that among all the 
configurations of equispaced concentric circular trajectories 
in Bn discussed in Section III-B3I the lowest path density is 
given by ^ achieved when A is equal to the upper limit of ^. 
A similar conclusion also holds for the union of interleaved 
spiral trajectories in Section III-B4I The lowest path density 
among all configurations of interleaved spiral trajectories in 
Bn discussed in Section III-B4I is also equal to ^ achieved 
when the parameter c is equal to the upper limit of — . Thus 



B. Optimality for where d > 3 

For fields in with > 3, we consider only unifom 
sets of the form ( l25l l. We also restrict ourselves to fields 
bandlimited to sets Vt that are compact convex subsets of 
and have a point of symmetry at the origin. As in Section 
III-Cl let Fn C Nn denote the collection of all uniform sets 
in K'^ that form Nyquist trajectory sets for Q,. From Theorem 
12.81 we know the necessary and sufficient conditions on the 
vectors {wi, U2, . . . , Wd} required for P e Fu. We now seek 
the solution to the problem 

min t{P) (34) 

where O C is a compact convex set with a point of 
symmetry at the origin. In this section we outline a procedure 
for solving the above problem. In our approach we relate this 
problem to the problem of designing optimal sampling lattices 
for static sampling in M''^^. 

Let P be a uniform set as defined in ( l25l l where p,„ is 
defined in ( l24l i with the vectors {ui, U2, ■ • • , vj} forming a 
basis for with (w^, Vd) = 5id for all i £ {1, . . . ,d}. Let U 
denote a d x d unitary matrix such that Uvd = e'', the unit 
vector along the d-th principal axis. Define 

tt^ {s e K'^-i : [/-I Q ^ G 17} (35) 

Also let Vi := Uvj and let Vi E M''^^ be the vector obtained 
from Vi by omitting its last component. We have the following 
result that relates the problem of designing optimal uniform 
sets to the well-studied problem of designing optimal sampling 
lattices. 

Proposition 3.4: Let C M'* be a compact convex set with 
a point of symmetry at the origin and let P be a uniform set in 
R'' as described above. Let U, Cl'-^ and Vi be as defined above. 
Then the uniform set P satisfies P G Afn if and only if the 
lattice of points defined by {^'iZi "i^A '■ ™ ^ Z'^^^} 

forms a sampling lattice for tl^ . Furthermore the path density 
£{P) is equal to the sampling density of b^ . 
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Proof: Consider the field / defined by 
f{r) = /(J7"V),r e W^. 
Clearly the Fourier transform of / is supported on the set 



{C/s ; s e O}. 



Now the problem of sampling the field / along the line pm {t) 
is equivalent to sampling the field / along the line Pm{t) 
defined by 



Pmit) 



d-1 

E 

1=1 



t e 



(36) 



where Vi = Uvi. Hence if P := {p,„ : m <^ TU^ ^} then it 
follows that 

P e AAn ^ P e M^-iu 

and that the path density of P is identical to that of P. 
We know from Theorem 12.81 that the necessary and sufficient 
condition for P e Mqu is given by 

\ ^ m,u^ i for all m e If^-^ \ {0} (37) 

4=1 

where Ui are defined as vectors in that satisfy {ui,Vj) ~ 
27TSij for 1 < i < (i — 1 and 1 < :/ < Now since Vd — e"^, 
it follows that {ui, e'^) = for i < d — 1. This fact together 
with the fact that fl is convex and symmetric about the origin 
implies that the condition (l37l i is equivalent to 



rf-i 
j=i 



^ , for all m € Z''~^ \ {0} 



(38) 



where Ui e ]R''-i is the vector obtained from Ui by omitting 
its last component. By construction it is clear that {ui,Vj) = 
{ui,Vj) = 2TTSij for I < i, j < d — 1. Thus it follows from 
lfT9l that condition ( |38] | is exactly the necessary and sufficient 
condition to ensure that forms a sampling lattice for Ct^ . 

We now consider the path density £{P). Since Vd = e'^ it 
is clear that the collection of points 

d-l 

m^Vi + tiid : m e Z"*, < e R} 

i=l 

remains unaltered if we replace Vi by 

■Si 


for 1 < i < d — \. Thus it follows via Lemma 12.91 that the 
path density t{P) satisfies 

i{P) = |det(G)ri 

where Gis ad— Ixd— 1 matrix with entries G,;.j = 
{vi,Vj) = {vi,Vj). Since the sampling density of b'^ is equal 
to |det(G)|~i (see HI) and €(P) = £{P) the result follows 
from Lemma 12.91 ■ 
As an immediate consequence of the above result we have 
the following corollary on optimal sampling trajectory sets 
from Ta- 



Corollary 3.4.1: Let C M'' denote a compact convex set 
with a point of symmetry at the origin. Among all possible 
choices of d x d unitary matrices let U be the one such 
that the set ft^ defined as in ( |35] | admits a sampling lattice 
with minimal sampling density in K'^^^. Also suppose that 
the vectors {wi,W2, ■ . ■ , Wd-i} C M''^^ generate an optimal 
sampling lattice for fields bandlimited to Cf' . Let P be a 
uniform set as defined in (l25T l where the vectors Vi are given 
by 



U' 







,1 < i < d — 1 and Vd — 



(39) 



Then P e and solved the optimization problem (|34] |. □ 

The problem of identifying sampling lattices with minimal 
density is well studied in literature (see, e.g., lfT9l . 1341 . 
|l33). Such results can be used in conjunction with the above 
corollary to design optimal uniform sets in M''. Below we 
present some examples of C M'' for which we can explicitly 
solve for the optimal uniform set. 

Example 3.1: Suppose il is the spherical ball C R''. In 
this case it is easy to see that d'-^ defined in (l35l l is a d — 1- 
dimensional spherical ball of radius p in R'^"^ for all choices 
of the unitary matrix U. Hence without loss of optimality 
we choose U in Corollary 13.4.11 to be the identity matrix. 
Now suppose that {wi,W2, . ■ ■ jWd-i} C M''"^ generate a 
sampling lattice with minimal sampling density for fields in 

Define 



"''-1 bandUmited to B"^- 







,1 < i < d- 1 



Vd 



Then it follows via Corollary 13.4.11 that the uniform set P 
defined in (l25l l with the above choices for vectors Vi achieves 
the minimum in ( l34l i. For d ~ 3 the vectors Vi can be chosen 

as 



Vl 





This follows from the fact that a sampling grid based on 
a hexagonal lattice is the optimal lattice (19\ for two- 
dimensional isotropic fields. This optimal configuration of a 
uniform set and the associated hexagonal sampling lattice is 
illustrated in Figure [TT] The exact choices of vectors Vi for 
fields bandlimited to spherical balls in R** for all d < 8 can 
be obtained from the above result using the results on optimal 
lattices for isotropic fields presented in fTPl Table C.I]. As 
mentioned in fT9ll . these results are based on results on closest 
packing of spheres in R'*. □ 

Example 3.2: Now suppose 51 is a cuboidal region in R'^ 
given hy 0. = {{x, y, z) : \x\ < p^, \y\ < py, \z\ < p^} where 

"^We are being imprecise here to keep the presentation simple. In reality the 
optimal orientation U can be found. But the optimal uniform set P can only be 
approached since an optimal sampling lattice for fj'^ can only approached. To 
be precise, one would have to consider a series of uniform sets P'^ satisfying 
lim£_>o (-{P"^) = ™fQs J"n ^(Q) akin to the statement of Theorem 13.21 The 
P mentioned here is the limit of P". 
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Fig. 11. Optimal configuration of a uniform set for sampling an isotropic 
field in . The cross-section forms a hexagonal sampling lattice as shown. 



Px < Py ^ Pz- We are interested in choosing the vectors 
vi,V2, V3 in a uniform set P of the form (IZSl l such that it 
solves (|34] |. In this case it can be seen that the mentioned in 
Corollary |3AT] admits a lattice with minimal sampling density 
when we have ~ e^. In this case U is the identity matrix 
and the set ^l^ = {{x, y) : \x\ < px, \v\ < Py}- Therefore, by 
the corollary it follows that (|34] | is solved by 



Vl = 



V2 = 




Such a trajectory set is illustrated in Figure |9] where Ai = — 



and A2 = — . 



□ 



In the following proposition we present a different opti- 
mality property of uniform sets. We show that uniform sets 
are optimal among all sets in Afn that contain all points in 
a sampling lattice for ft. This result is the R'* analogue of 
Proposition 13.31 for K^. 

Proposition 3.5: Let C M'^ denote a compact convex set 
with a point of symmetry at the origin, and Q G Afn be any 
trajectory set that visits all points in a sampling lattice for fl. 
Then there exists a uniform set P e Afn such that £{P) < 

my 

Sketch of proof: This follows from Lemma 13.11 and 
Lemma l272l via the same argument used in proving Proposition 
[33] ■ 
The above result implies that the optimal uniform set given 
by Corollary 13.4.11 is optimal among all trajectory sets that 
visit all points in a sampling lattice for fi. However, we 
have not solved ( [29l ). Nevertheless, as we mentioned in the 
case of trajectories in M^, it may be possible to extend 
Landau's necessary conditions ||3T1 and the results on optimal 
configurations of points 1321 , ||33l Cor 4.5] for sampling in 
R'' to Nyquist trajectory sets in R"* and thus obtain a solution 
to 

IV. Generalization to higher dimensional 

SAMPLING MANIFOLDS 

In Sections and |III] we studied the problem of sampling 
bandlimited fields in R'^ on trajectories. Trajectories can 
be regarded as one-dimensional manifolds in R''. We now 
consider a generalization of this problem in which we replace 
trajectories with general K-dimensional manifolds in R'' for 
K < d. Such a setting is relevant in applications where 



the process of measuring a d-dimensional field involves the 
process of recording the field values on several K-dimensional 
manifolds at a fine resolution. For example consider the prob- 
lem of sampling and reconstructing a 3-dimensional object 
from its views along several 2-dimensional sections. Such 
a sampling and reconstruction scheme is often employed in 
applications like Transmission Electron Microscopy (TEM) 
Is] and Magnetic Resonance Imaging (MRI) f9l where one 
tries to recreate the 3-dimensional structure of an object after 
imaging or scanning the surfaces of several cross-sections of 
the object. In such cases, it may be of interest to identify the 
optimal orientations along which the images should be taken 
so as to minimize the total area of the images taken while still 
being able to reconstruct the 3-dimensional structure of the 
object. 

We use the following terminology. A n-manifold in 
will mean a topological manifold M embedded in R'' such 
that every point in 1\I has a neighborhood in M that is 
homeomorphic to R". E.g., a 1-manifold in R'^ is a curve 
in R'' and a 2-manifold in is a surface in R'^. Generalizing 
the concept of a trajectory set, we define a K-manifold set 



as a countable collection of K-manifolds in 



of the 



form P 



G 1} where pi denotes a K-manifold 



in R and I C Z. We also define a generalization of the 
path density metric to manifolds. If one assumes that it is 
inexpensive to increase the sampling density on the sampling 
manifolds, the metric of interest is the total volume of the all 
the sampling manifolds per unit spatial volume. We call this 
metric the manifold density of the manifold set. For a set P 
of K-manifolds in R"^, the manifold density is defined by 



PniP) ■= lim sup 



sup 



Void (a) 



(40) 



where V^(a, a;) represents the total K-dimensional volume of 
the portions of the manifolds from P that are located within 
the ball B^{x) and Volci(a) represents the volume of the ball. 
For example, when k = 1 the manifold density is exactly the 
same as the path density defined in (|5]l, and when k = 2 
and d — 3, the manifold density P2{P) represents the total 
area of the sampling surfaces per unit volume. Similarly, we 
generalize Nyquist trajectory sets to manifolds. 

Definition 4.1: A K-manifold set in R'' of the form P = 
{pi : i G 1} is called a Nyquist K-manifold set for il C R'' if 
it satisfies the following conditions: 

(Ml) [Nyquist] There exists a uniformly discrete 
collection A of points on the manifolds in P 
such that A admits perfect reconstruction of 
fields in Bn, i-C-, A C {r : r € pi, i € 1} and 

A e An- 

When the value of k, is clear from context, we just say Nyquist 
manifold set. We use a special notation for collections of 
Nyquist K-manifold sets: 

Definition 4.2: For fi C R'^ we define Af^ as the collection 
of all Nyquist K-manifold sets for Q,. 

The examples of Nyquist trajectory sets in R'' that we consid- 
ered in Section nil are Nyquist 1-manifold sets, or equivalently, 
elements of Af^ for k — 1. We now seek generalizations of 
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these examples to k > 1. The other problem of interest is 
a generalization of ( |29] ). However, we note that the problem 
minpgAA" Mk(^) is in general ill -posed under the current defi- 
nition of AfJ since it is possible to construct a set of manifolds 
with vanishing manifold density, e.g., the degenerate case in 
which each manifold pi is confined to a small neighborhood 
of some point in a sampling lattice for il. This problem can 
be made more meaningful if further restrictions are placed 
on the elements of Af^ by imposing a condition analogous 
to condition (C|2) in the definition of J\fn in Section III-AI 
However, we do not consider such a definition in this paper 
to keep the presentation simple. Instead, restricting ourselves 
to hyperplanes, we establish a generalization of Theorem 13.21 
for sets of hyperplanes, i.e., affine manifolds of dimension 
K = d — 1. 

For some h G with \\h\\ = 1, let i7 denote the 
hyperplane H := {r G M."^ : {r,h) =0}. Generalizing uniform 
sets of straight line trajectories we define a uniform hyperplane 
set in K'^ to be a collection of equispaced parallel hyperplanes 
of the form P ~ {pi : i S Z} where 

Pi = w + iAh + H 

is a shifted versioijf] of the hyperplane H. Thus P is a 
collection of hyperplanes parallel to H located at intervals 
of A units. The manifold density of the uniform hyperplane 
set P is given in the following result. 

Lemma 4.1: The manifold density of the uniform hyper- 
plane set P defined above is given by 

where A is the spacing between adjacent parallel hyperplanes. 

Sketch of proof: This lemma can be proved in exactly 
the same way as Lemma [Z!2] bv considering d— 1-dimensional 
rectangular regions located within a d-dimensional spherical 
ball. We skip the details. ■ 
We first obtain generalizations of Proposition 12.31 and Theo- 
rem 123] to unions of uniform hyperplane sets. Given finite sets 
of vectors {wi\f^^, unit vectors {hi}f^^, and scalars { A^jfl]^, 
let Pi denote the uniform hyperplane set 



P^ = {P^.J : J e Z} 



where 



jA,h, + H, 



(41) 



(42) 



represents a shifted version of hyperplane i/; := {r e R'^ : 
(r, hi) = 0}. The union of the N uniform hyperplane sets 
P, : 1 < i < is defined as 



N 



P := 



(43) 



We are interested in finding the conditions on P so that P e 
A/"^"^. Since we assume that the field can be sampled at a 
fine resolution on each of the sampling manifolds Pi, we can 
construct a sufficiently fine lattice in R'' composed of points 
from Pi such that every point on Pi is arbitrarily close to 

^We avoid using tlie notational convention of (5) liere wliicli we reserve 
for sets in tlie frequency space. 



some point from the lattice. For such a lattice, it follows via 
classical sampling theory ||T9l that the spectral repetitions in 
the sampled spectrum obtained from the samples taken on such 
a lattice are spaced arbitrarily far apart in all directions except 
one. In this case if is a compact set, the sampled spectrum 
from observations on Pi satisfies: 



c^p{\{jui,Wi))F{uj + jui),uj e n (44) 



where Ui := This is the analogue to ( fTSl l and is derived 



explicitly in Lemma IA.41 in Appendix |K 
Let Q cW^ denote the set of points 




: k, e {0, 1},1 < i < 



Then we have the following results that generalize Proposition 
12.31 and Theorem 12.41 The following proposition gives a 
necessary condition that must be satisfied by any union of 
uniform hyperplane sets that forms a Nyquist manifold set for 

n. 

Proposition 4.2: Let C K.'* be a compact convex set and 
let P denote a union of uniform hyperplane sets of the form 

dllll. Suppose Q C 0(s) for some s e R'^. Then P ^ A^^^^ 

□ 

As in the case of Proposition 12.31 the above proposition can be 
proved by constructing a sinusoidal field that vanishes on all 
the hyperplanes in P, and has a Fourier transform supported 

o 

on Q(— s) C il. We provide a proof in Appendix |J] In other 
words, this proposition implies that every P e 
necessarily satisfy 



must 



Q ^ 0(s) for all s e R''. 

This necessary condition must be satisfied by all uniform hy- 
perplane sets P that form Nyquist manifold sets for fl. When 
the vectors hi are linearly independent (which automatically 
forces the constraint N < d), it can be shown that the tightest 
necessary condition given by Proposition l4.2l is also sufficient. 
In that scenario, the following theorem provides necessary and 
sufficient conditions to ensure that P forms a Nyquist manifold 
set for O. 

Theorem 4.3: Let 17 C R'' be a compact convex set. Let 
P denote the union of uniform hyperplane sets defined in 
(|43] |. Suppose that the vectors {hi, /12, . . . , /iat} are linearly 
independent. Then we have the following generalization of 
Theorem 12.41 



P e 7Vj2~^ if Q ^ n{s) for all s e R"*, (45) 
and moreover, 

P i M'l^^ if Q C h{s) for some s G R'^. (46) 

□ 

We provide a proof in Appendix |K] The proof of the first 
result (|45l l is very similar to that of the analogous result 
(fT9] l for R^. We also note that the second result ( |46] l is a 
restatement of Proposition 14.21 As we saw in the discussion 
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following Theorem 12.41 these results can also be interpreted 
as the conditions for sampling on a union of shifted lattices 
1241 under special case in which each of the lattices is finely 
sampled in all directions except one. We do not repeat the 
discussion here since the argument is very similar to that for 
uniform sets in M^. 

Below, we illustrate the theorem using a simple example of 
sampling a field in R^. In such a case, the hyperplanes in M'' 
are just planes. 

Example 4.1: Suppose O C M"^ is a spherical ball of radius 
p centered at the origin. For A,; > let Sj-s = (-^,0,0), 
Sy^i = (0, f^,0), and s^,, = (0,0, f^-). We consider three 
examples of unions of uniform hyperplane sets. 

1) Uniform set of planes Consider a uniform set of planes 
parallel to the yz plane spaced at intervals of Ai units 
given by P = {pi : i £ Z} where 

p, = {{iAi,y,z) :y,ze R}. 



Such a configuration is illustrated in Figure [T2( a) | In this 
case, we have spectral repetitions along uJx axis. Hence 
it follows from Theorem |43] that the maximum value of 
Ai to ensure that P £ J^q~^ is A| = ^. It is easily 
verified via elementary geometry that Ai < A^ is the 
condition to ensure that the intersection of two spheres, 

n il.{sx,i) is empty. 

It follows that this is exactly the condition to ensure that 
there is no aliasing in ( l44l i. 
2) Union of two orthogonal uniform sets of planes Consider 
a union of two uniform sets of planes one parallel to 
the yz plane spaced at intervals of A2 units and another 
parallel to the xy plane spaced at intervals of A2 units 
given by P = Uf^^Pi where Pi = {pi,j : j € Z} where 

Pi J = {0'^2, y,z):y,z£ R} and 
P2j = jAz) : a;,y e R}. 



Such a configuration is illustrated in Figure |12(b)| In 
this case, we obtain spectral repetitions along the lj^ 
axis from Pi and along uj^ axis from P2. Hence it 
follows from Theorem 14.31 that the maximum value of 



A2 to ensure that P e M^f^ is Aj = It is easily 
verified that A2 < A2 is the condition to ensure that 
the intersection of three spheres, 

n ^{sx,2) n ^{sz^2) is empty 

In other words, for each w G there is at least one 
set among 0,{sx^2) and 51(5^.2) that does not contain 
Lo. It follows via (|44] | that this is exactly the condition 
required to ensure that for every w G fi, there is at least 
one i £ {1,2} such that the sampled spectrum P](i^) 
from the samples on Pi matches the true spectrum P(w). 
3) Union of three mutually orthogonal uniform sets of 
planes Consider a union of three uniform sets of planes 
one parallel to the yz plane, another parallel to the xy 
plane and a third parallel to the xz plane, with each 



spaced at intervals of A3 units given by P = uf^^P^ 
where Pi = {pi.j : j E Z} where 

pij = {(jA3,?/,z) : e R} 
P2j = {(a;,2/, jAa) : a;,y e M} and 
P3.3 = {{x,jA3,z):x,zeR}. 



Such a configuration is illustrated in Figure [T2(c)| In this 
case, we obtain spectral repetitions along the lUx axis 
from Pi, along lj^ axis from P2 and along ojy from P3. 
Hence it follows from the theorem that the maximum 
value of A3 to ensure that P e A/"^"^ is A3 = 
It is easily verified that A3 < A3 is the condition to 
ensure that the intersection of four spheres, 

il n n{sx,3) n fi(s.y^3) n il(s2,3) is empty. 

In other words, for each w G there is at least one 
set among il{sx,3), ^{sz.s), and ^(5^,3) that does not 
contain oj. It follows via ( l44b that this is exactly the 
condition required to ensure that for every lj £ fl, 
there is at least one i £ {1, 2, 3} such that the sampled 
spectrum Fs{uj) = F{uj), the true spectrum. 

□ 

We now present a generalization of the optimality result of 
Theorem [32] We denote by £^ the collection of manifold sets 
in M^^^ that are unions of uniform hyperplane sets. We have 
the following optimality result for a uniform hyperplane set. 

Theorem 4.4: For any nonempty compact convex set Vl C 
R'', let u be the minimizer in ( |32] |. and for e > let P*^ denote 
a uniform hyperplane set given by P*^ = {pj : j S Z} where 



(- 



27r 



u 



with Hu = {r £W^ : {r,u) = 0}. Then P' £ £^ for all 
e > and is optimal in manifold density as e — > 0, i.e.. 

Win) 



lim /i«(P') 



inf 



27r 



□ 



This result is established by identifying the union of uniform 
hyperplane sets with minimal manifold density that satisfies 
the conditions of Theorem 14.31 We provide a proof in Ap- 
pendix |L] In short. Theorem 14.41 establishes the optimality of 
a uniform set P from all trajectory sets in £q. In particular, we 
have identified a sequence of uniform sets of hyperplanes in- 
dexed by e with manifold densities converging to the infimum 
manifold density from £^ as e tends to zero. The limiting 
uniform set corresponding to critical sampling is obtained by 
plugging in e = 0. This critically sampled uniform hyperplane 
set is a set of parallel hyperplanes oriented parallel to u and 
spaced apart. 

For fields in M'', this result implies that for sampling a 3- 
dimensional field it is more efficient to sample along one set of 
equispaced parallel planes rather than along several different 
sets of equispaced parallel planes. In particular, in Example 
14. 1[ this result implies that the first configuration of planes 
has a lower manifold density than the other two examples. 
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(a) Uniform set (b) Two orthogonal uniform sets 

Fig. 12. Three configurations of samphng on planes discussed in Example 14.11 
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(c) Three orthogonal uniform sets 



This conclusion can be verified by calculating the optimal 
manifold densities in the three cases using Lemma 14.11 The 
minimum manifold densities in the first, second and third cases 



are respectively, 



x/2p 



and 



^/3p 



V. Reconstruction schemes 

We now consider schemes for reconstructing bandlimited 
fields using measurements of the field taken along the various 
sampling trajectories and manifolds proposed above. The 
reconstruction schemes we propose can be broadly classified 
into two types. For trajectory sets and manifold sets composed 
of affine sets such as lines and hyperplanes, we use results 
from sampling on lattices and unions of lattices. For non-affine 
trajectories and manifolds, we use schemes for non-uniform 
sampling based on the Beurling frame theorem ifTSl . 

We first consider sampling trajectories. We know from 
Lemma 12.11 that uniform sets of the form (|7]i for fields in 
and (ESll for fields in R'' satisfy condition (iffl. Hence 
the field values at all points on these straight line trajectories 
can be reconstructed from samples taken at uniform spatial 
intervals. The maximum spacing allowed between adjacent 
samples on the lines can be calculated from Lemma 12.11 
Suppose that the sensors on each of the parallel lines take 
samples at uniform intervals and that the sample locations on 
the various parallel lines are aligned with each other. In such 
a case the collection of all samples obtained on all these lines 
effectively corresponds to a collection of samples of the field 
taken over a periodic lattice of points. Hence, for uniform 
sets of trajectories, any reconstruction algorithm used for 
reconstructing bandlimited fields on lattices is sufficient (see, 
e.g., |fT9l ). For unions of uniform sets like in (fTST i the resulting 
set of points forms a union of sampling lattices. In this case 
reconstruction schemes for sampling on unions of lattices 
are applicable (see, e.g., |[24l). Alternatively, more general 
schemes for non-uniform sampling such as the projection on 
convex sets (POCS) algorithm 1361 can be used. 

Reconstruction schemes for the circular trajectories consid- 
ered in Section III-B3 1 are provided in |[l4l and ITSl . However, 
since these trajectory sets do not satisfy condition (KH), these 
reconstruction schemes require the exact field values at all 



points on the circular trajectories. For an approximate recon- 
struction from samples taken by sensors moving at constant 
angular velocity one can use the approximation provided in 
(l22l l in conjunction with the reconstruction scheme of ifTsl . 
In addition, if one is interested in reconstructing the field 
accurately only over a finite number of lines through the 
origin, then it is sufficient to sample the field at all points 
of intersection between these lines and the circles and use the 
reconstruction scheme of ifTsl . Alternatively, as we showed 
in Section III-B3I it is possible to reconstruct the field from a 
uniformly discrete collection of points on the concentric circles 
chosen to satisfy the conditions of Beurling's covering theo- 
rem. The reconstruction algorithm is based on a frame analysis 
and is described in 1181 . The same reconstruction algorithm is 
also applicable for the union of spiral trajectories of Section 
III-B4I Similarly, even for arbitrary non-affine trajectories like 
the one shown in Figure |l(b)| it is possible to identify a set 
of sampling points and a corresponding exact reconstruction 
scheme via the Beurling theorem. 

Now consider the case of sampling on uniform hyperplane 
sets. We can make an argument analogous to that of sampling 
on uniform sets of trajectories. The restriction of a bandlimited 
field to the hyperplane is also bandlimited. Like in (|42] | any 
hyperplane H in can be expressed in the form 

i7 = {r + 6 : (r, n) = 0, r e M''} 

for some vectors n, e W^. Now if we let A denote a d x 
d — 1 matrix with linearly independent columns that are also 
orthogonal to n, then the hyperplane H can be expressed as 
H = {Ax + b : X £ W^^^}. In the following lemma we 
identify the bandwidth of the field / restricted to H. 

Lemma 5.1: Let {/(?') : r G M'^} denote a field bandlimited 



to some set ft C 



Consider any hyperplane H C 



expressed in the form H = {Ax + b 



X e 



^} where 



Aisa.dxd— 1 matrix. Then the restricted field g defined by 



gix) = fiAx + b), xe 



is bandlimited to fie := {A'^ 



en} c 



pd-l 



Proof: Using F and G to denote the Fourier transforms 



of / and g respectively, we have 

g{x) — I F(tt') cxp(i(tt', Ax + 6))(iw 

whence 

G(s) = / F{Lu)exp{\{{uj,Ax + b)-{s,x)))dxdio 

= I F(w)exp(i(w,fo))(27r)''~^(5(A^a;-s)da; 
= for s ^ r^G- 

■ 

The conclusion of the lemma above shows that the restriction 
of a bandlimited field to a hyperplane is also bandlimited. 
Hence the field values on all points on the hyperplane can be 
reconstructed exactly using only the field values measured on 
a sufficiently fine lattice of points on the hyperplane. Thus a 
sampling scheme using a uniform hyperplane set of the form 
in W\\ can be practically implemented by using measurements 
on a d — 1 dimensional lattice of points for each of the hy- 
perplanes. If the sampling lattices on each parallel hyperplane 
are aligned with each other then the collection of all sampling 
points from all the hyperplanes forms a lattice in M"^. Thus 
reconstruction schemes for reconstructing bandlimited fields 
on lattices can be used in these cases (see, e.g., |fT9l , 1361 ). 
Similarly, if we consider sampling manifold sets like in (l43T l 
composed of unions of uniform hyperplane sets, reconstruction 
schemes for sampling on unions of lattices are applicable (see, 
e-g-j 1241 '). Alternatively, the POCS algorithm [I36| can also be 
used for non-uniform sampling. 

VI. Conclusion and future work 

In this paper we have introduced the problem of sampling 
bandlimited fields in on trajectories and other lower 
dimensional manifolds. We have presented various examples 
of Nyquist trajectory sets and Nyquist manifold sets that 
admit perfect reconstruction of bandlimited fields. We obtained 
necessary and sufficient conditions on a union of uniform 
hyperplane sets to be able to reconstruct bandlimited fields. 
We also introduced the path density and manifold density 
metrics and illustrated some optimality properties of select 
trajectory sets and manifold sets in terms of these metrics. In 
particular, we established the fact that a uniform hyperplane 
set achieves the minimum manifold density among unions of 
uniform hyperplane sets that admit perfect reconstruction of 
bandlimited fields. Some of our main results are summarized 
in Table U Besides these results, we also have Propositions 
13.31 and 13.51 which show that uniform sets in and R'^ are 
optimal from among all Nyquist trajectory sets that visit all 
points in a sampling lattice for a given Vl. In addition, we have 
considered non-affine trajectories such as concentric circles 
and interleaved spirals and discussed some known results 
and some new results on sufficient conditions for perfect 
reconstruction. 

This paper opens numerous avenues for future work in terms 
of extensions and generalizations. We have not completely 
solved the general problem that we posed in (|29] | which seeks 
the trajectory set in Mn with minimal path density. The best 
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Lines 


Hyperplanes in 


Unions 

sets in M 


of uniform 

2 


Uniform sets in R"* 


Unions of uniform hy- 
perplane sets 


Prop 12.3 


(Nec) 


Tlim|Z8j(Nec & Suf) 


Prop 14.2 


(Nec) 


TlimlM 


(Nec & Suf) 


ProplMKNec & Suf) 


Thml4.3 


(Nec & Suf) 


Tliml3.2 


(Opt) 


CorB.4.1l(Opt) 


Thi-nl4.4 


(Opt) 



TABLE I 

Summary of some results. Each result uses the one above it. 

"NEC" stands for necessary conditions and "Suf" for 
sufficient conditions for perfect reconstruction and "Opt" 
stands for optimality conditions under the sampling scheme 
specified by the column. 



results we have so far are Theorem 13.21 Corollary 13.4.11 
and Proposition 13.31 which establish different sub-optimality 
properties of uniform sets. We also made the interesting 
observation that for isotropic fields bandlimited to a circular 
disc Vl of radius p in M?, the minimum path density for 
uniform sets in Mn matches the minimum path density for 
the sets of concentric circles that we considered as well as 
the minimum path density for the set of interleaved spirals we 
considered. This minimum value is given by It remains to 
be seen whether this can be bettered by employing some other 
trajectory sets, perhaps involving a non-uniform collection 
of sampling trajectories. As mentioned earlier, a potential 
approach would be to try to generalize Landau's lower bound 
on the minimum sampling density required for sampling 
bandlimited fields on points |f3T| (see also |!20! Corollary 1]) to 
analogous bounds on the minimum path density and minimum 
manifold density required for sampling on trajectories and 
manifolds. It may then be possible to extend known results 
on sampling configurations with minimal sampling density 
l32l . 1331 Cor. 4.5] to Nyquist trajectory sets to solve ( |29l ). 
Similar extensions to the problem of designing manifolds 
with minimal manifold density are also of interest. However, 
before proceeding we need to first generalize condition (C® 
to manifolds in such a way that the class of allowed sampling 
manifolds can be restricted to a meaningful set that does not 
admit trivial solutions composed of degenerate manifolds. 

Many results in this paper can potentially be extended to 
more general settings. For instance, most of our results are 
under the assumption that the field is bandlimited to a convex 
subset of R''. Ideas from the proofs of these results can 
be used to extend these results to more general classes of 
fields, e.g., bandpass fields bandlimited to non-convex regions, 
and non-bandlimited fields that form shift-invariant spaces. 
The results on sampling on circles and spirals in R^ can 
be extended to analogous results on sampling on concentric 
spherical shells, concentric cylindrical shells, and concentric 
helixes in R'' for d> 2. The Beurling theorem can be applied 
to design sampling and reconstruction schemes for such non- 
affine manifolds. 

An important practical aspect that we have ignored in this 
work is ambient noise which affects the measurement process. 
In the presence of noise, perfect reconstruction is not possible. 
In such a case, it would be of interest to study the tradeoff 
between path (manifold) density and SNR in the reconstructed 
field and to identify optimal sampling trajectories (manifolds) 
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that optimizes the tradeoff at a given noise level. This tradeoff 
is also relevant in sampling non-bandlimited parametric fields 
such as diffusion fields using mobile sensors. 

It is also of interest to extend our analysis to the design 
of fc-space trajectories for MRI. The bandlimited field model 
we use in this paper is analogous to using a spatially limited 
field model in MRI. However, in MRI literature one typically 
assumes that in addition to the object being spatially limited, 
the object has more energy in the lower frequencies in /c-space. 
It would be interesting to generalize our approach to such a 
setting wherein one makes assumptions on the field both in 
the spatial and frequency domains. 
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A. List of symbols 
/ : C 
F ■.R'^^C 

\M 

{v,w) 
S{.) 

Uniformly 
discrete set 
Ba 

An 

Bi{x) 



Appendix 
and definitions 
Field of interest 

Fourier transform of the field of interest 
defined in ([T]i 

Euclidean norm of the vector v 
Euclidean inner product between vectors v 
and w 

Dirac delta function 

Kronecker delta function evaluated at (i, j) 

Any set A C M"^ satisfying inf{||x — yj| : 

x, ?/ e A, x 7^ y} > 

Collection of fields bandlimited to a given 

set Vl as defined in ^ 

Collection of all uniformly discrete sets 



A c 



with the property that any field 



Uniform set 
in m? 

Uniform set 

in R'^ 



/ e So can be reconstructed exactly from 
its values on A 

Set VL shifted by the vector s as defined in 
© 

Path density of trajectory set P defined in 
© 

The d-dimensional spherical ball of radius 
a centered at the origin in W'-; := {y E 
K'^ : ||y|| < a} 

The d-dimensional spherical ball of radius 
a centered at x e W^: B^{x) := {y G K'^ : 
\\y-x\\ < a} 

Collection of all Nyquist trajectory sets for 
il, i.e., trajectory sets that satisfy condi- 
tions (ClD and (dD 

Trajectory set composed of equispaced 
parallel lines as defined in © 
Trajectory set composed of a periodic con- 
figuration of parallel lines of the form in 



£n 



Collection of all uniform sets that are in 
A/n 

Collection of all trajectory sets in Mq, 
that can be expressed as a finite union of 
uniform sets 

Width of a convex set as defined in ( l32l i 
Manifold density of a set P of K-manifolds 
as defined in ( l40l i 

Collection of all Nyquist K-manifold sets 
for f2, i.e., manifold sets satisfying condi- 
tion (IV© 

Collection of all manifold sets in M^^^ 
that can be expressed as a finite union of 
uniform hyperplane sets 



B. Proof of Lemma \2.2\ 

Let Pi be the straight line trajectory defined in ©. Let £i 
denote the length of the portion of the line pi located within the 



disc Bl{x). Let 



denote the range of indices i such that 



ii > 0. We are interested in the total length of the portions 
of the lines within the disc B^{x). This total length is as 
illustrated in Figure |13(a)| when the lines are horizontal. Let 
i = arg maXj £i . Now we can approximate the area of the disc 
above and below by rectangular regions with sides £i x A to 
obtain the inequality: 



£,A - £jA < Tra^ < ^ £,A 



' 2a A. (47) 



The areas representing the lower and upper bounds are 
illustrated respectively in Figures |13(b)| and |13(c)| Now 

T>^{a,x) = ^i- Hence we have 



VP{a,x)-£. 



< 



< 



VP{a,x) +2a 



and (fT2l i follows by taking hmits. It is also immediate that line 
segments of length £j for i G can all be connected by a 

continuous curve by connecting portions the ends of adjacent 
segments via a part of the circumference. The length of such 
a curve is no more than 2?^ (a, x) + ira and hence P satisfies 
condition (C©. 



C. Proof of Proposition 12.31 

o 

By the given condition it follows that Q{—s) C fi. Now 
consider the field 

/e(r) := ccxp(-i(s,r))nsin ( i^^Hl,^^ ,r G 



The Fourier transform of the field is supported on the set 
Q{—s) in the Fourier domain. Although the fields are not 
in L^(C) they can be approximated by functions G L^(C) 
defined by 

2 

9c{r) = /c(r) J|sinc(— ) 
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(a) Total length of interest (b) Area representing lower bound in 




(c) Area representing upper bound in 

ED 

Fig. 13. Computing the path density of a uniform set 



where sinc(a;) 



sin(7ra:) 



and Ti denotes the i-th component 



of vector r. The Fourier transform of gc is supported on points 
close to Q(— s). In particular, we have 

G,H=Oifw^UteQ(_,)i?,(0 

where i?e is the square [— e, e]^ and i?c(t) is the set i?e shifted 
by t following our convention from ([3]l. 

Since e can be made arbitrarily small, it follows from the 

o 

fact that Q(— s) C that there exists an e such that Gc is 

o 

supported on a subset of the set il in the Fourier domain. It 
can also be verified that for all c G C, the functions gc vanish 
on the straight line trajectories in P. Hence, all fields gc are 
bandlimited to O but they cannot be reconstructed uniquely 
based on their values on the lines in P. Thus it follows that 
P does not satisfy condition (Cjl) and hence P ^ A/q. 



D. Proof of Theorem \2.4\ 

We prove only ( fT9] l since result of ( |20] i clearly follows from 
Proposition |23] We need the following definition and lemma. 
The lemma is made more general than what is required for 
this proof because the general result is required later 

For any set S C Z^, we say that S is lattice-convex if the 
following condition holds: if a, 6 S 5 then 5a + (l — 
for all (5 e (0, 1) such that 5a +{1- 5)h S Z^. For the lattice 
we define the unit-cell as the subset of points 



U 



N 



' N 



e {0, 1} for alH ^ C Z 



N 



where Xi denotes the i-th component of x, and e* denotes 
the unit-vector along the i-th principal axis. A translate of the 
unit-cell is any shifted version of the expressed in the 
form {x + z : X & U^} for some z G Z^. 

Lemma A.l: Let N C be a compact lattice-convex set. 
Associated with each element n e N is a value v{n). These 
values satisfy the linear relations 



9n 



E 



^v{m) l<i <N 



(48) 



where Ti^m G C \ {0} and and r^.m are known, and 

M„,i := {m G N : m — n = fee*, for some k £ Z}. 

Suppose that these equations are consistent, i.e., they admit 
at least one solution. Further suppose that N does not contain 
01 its translates. 
Then the values {v{n) : n G N} can be uniquely identified 
from the relations in ( |48] |. 

Proof: Since the equations in ( |48] | are known to be 
consistent, we just have to show that these equations admit 
a unique solution. This is equivalent to showing that the 
homogenous system obtained from those equations with the 
(7^ set to zero admits a unique solution. In other words we 
need to establish the following claim: 



Claim: Under the conditions of the statement of Lemma lA.ll 
the equations 



= 



E 

rneM„ 



iv{m) l<i < N 



(49) 



admit a unique solution given by v{n) ~ for all n e N. 

We prove the claim using induction on N. Suppose = 1. 
In that case it is easy to see that if N does not contain 
or its translates, N is a singleton and hence the claim holds 
trivially. Now suppose that the claim holds for N = — 1. 
Now let N = N'^. Suppose N does not contain or its 
translates. Define 



c{njvo : n = (tii, 712, . . . , n^o)^ G N} 



and 



Af(N) := |{njvo : n = {ni,n2, ■ ■ ■ .tlmo] 

where n^o denotes the A^^-th component of vector n. Let 

:= {n = (ni, 7^2, . . . , rif^o)^ G N : ripfo = V}. 

Now let C be defined as 

N2 {ne : n-e^" G N}. 

Clearly, by construction, for all m G \ N^, we have 
Mm.Afo = {"^}- Since r^ro „j ^ this implies that 



v{m) 



9: 



N" 



for ™ G \ 



and hence all values {v{m) : m G \ N^} are necessarily 
equal to 0. 
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It is easy to see that there is a set C such that one 

can define a bijection between elements n e and n e 
as 

n (ni,n2, . . . ,ri^,o_i)^ 

and n = (rii, 712, ■ • ■ , ri7vo-i> J^)^ = «2, • ■ • ,_nAro_i, J^)^. 
It is also easy to see from this definition that is lattice- 
convex, because N was itself lattice-convex. Furthermore, 

does not contain or its translates, since N does 

not contain or its translates. To each element fi G 

we associate a value := v{n) where n = 

{ni,n2, ■ ■ ■ , ri^o_i, 17)^. It is clear by the definition of and 
■w{.) that the equations in ( |49] l corresponding to ?i € and 
i € {1, 2, . . . , N'^ — 1} involve only the values {v{m) : m £ 
N^} and hence can be rewritten with v{m) replaced by w{rh) 
where rfi is the element of corresponding to m. Thus it 
follows via the induction assumption for = 7V° — 1 applied 
to A^^ that the values {w{h) : n e N^}, or equivalently the 
values {v{m) : m £ N^}, are necessarily equal to 0. Thus we 
have shown that v{m) = for all m £ N^. 

In order to complete the proof of the claim, we need to 
show that v{m) = for all jti e N \ N^. This follows by 
a straightforward induction on 7\/(N). Since = N for 
M(N) = 1, we have established that the claim is true for 
all lattice-convex sets N C Z^" such that A/(N) = 1. Now 
with N still fixed at iV = iV" assume that the statement 
of the claim is true for all lattice-convex sets N C 
with M(N) = — 1. Now consider any lattice-convex set 
N C with i\/(N) = By the argument presented in the 
previous paragraph, we have v{m) = for all to e where 

:= {n € N : n^qa ~ 77} as before. It is straightforward 
to see that N \ is a lattice-convex subset of with 
M(N \ N^) = — 1. Furthermore, since we have v{m) = 
for all m G the variables {v{rn) : m £ N^} can be 
removed from the linear relations ( |49] l to obtain linear relations 
involving only the variables {v{m) : m £ N \ N^}. By 
the induction assumption on M ~ — 1, it follows that 
v{m) = must hold for all m e N \ N^. This completes the 
induction on A/(N) and thus the induction on N. Therefore 
this completes the proof of the claim and hence the lemma. 

■ 

We now proceed to the proof of (fT9l l. Suppose 

Q ^ n{s), for all s £ M^. (50) 

Let F denote the two-dimensional Fourier transform of any 
field / bandlimited to ft. Let uj £ flhe arbitrary. Consider the 
following set 

N -.^ {n£Z^ ■.uj + Un£n} (51) 

where U is a 2 x N matrix with i-th column given by Ui 
defined in ( fT4b . It is easy to see via the convexity of ft that N 
is a lattice-convex subset of . Furthermore, condition ( fSOl l 
implies that translates of the unit-cell are not contained 
within N. To each element n of N we associate the value 
v{n) = F{uj + Un). Now the expressions in dTSI l for the 
spectrum F] of the samples taken on points on the parallel 



lines in Pi can be rewritten as 

exp(i(a;,w,))i^;(w) 

= ^exp(i(a; + jMi,Wi))i^(a; + e 17, (52) 

for all i £ {1, 2, . . . , iV}. Evaluating these equations at the 
points {uj + Un : n £ N} we obtain consistent linear equations 
in v{.) of the form ( |48] | with = cxp(i(cj + [/n, u'i))i^*(w + 
Un). Applying Lemma lA. II to the set N defined in ( fSTl i we 
conclude that for each lo £ Q. the values of the Fourier 
transform {F{ijJ -\- Un) : n £ U} can be decoded using the 
values {Fl{iu + Un) : n G N,l < i < N}. Since bj £ 
was arbitrary this proves that F{lo) can be recovered for all 
Lo £0,. Thus we have verified that P satisfies condition (Cjl) 
whenever (ISOl l holds. 

We now verify that P satisfies condition (CS). For each i, 
it is clear that for any x G the portions of the lines in the 
uniform set Pi within (x) can be connected by a continuous 
path containing just these portions as well as portions of the 
circumference of the circle. Thus the length of such a curve 
is just T)^^ (a, x) + 0(a). Hence N such curves corresponding 
to each of the {Pi]i can be joined into a single curve with 
length E!Ii 25^' (a, 2:) + 0(a) ^ V^{a,x) + o{a^). Thus 
condition (Cl2) in the definition of A/n also holds and hence 
P £ Afn- This completes the proof of (fT9] l. 



E. Proof of Lemma \2.6\ 
We have 

If. 

s{t) ~ — I F{ijj) cx^{\{aUxCO&(vt) + aUyS,in{h't)))dLjj 

where F is the Fourier transform of /. The Fourier transform 
of s is given by 

S{0 = ^ I F{Lu)K{Lu^,u:y,Odio (53) 
27r Jr2 

where 

= J exp(i(aa;a; cos(i^t) + aujy sin(i't) — £^t))dt 

= J cxp{\{Asm{iyt + B)-^t))dt (54) 

/+00 
JniA) cxp(i(n(i^t + B)- ^t))dt 
n— — oo 

+ 00 

= 27r ^ Jn{A)S{^ - ni')cxp{\nB) 

n— — 00 

where J„(.) denotes the Bessel function of order n and in 
( |54| ) we use A a||w|| and B := tan"^ Using the 
approximation that Jn{x) w for |n| > |.t| + 1 we get the 
desired result from ( l53T l via the bandlimitedness of F. 
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F. Proof of Lemma \2.7\ 

From the geometry of the problem it is clear that the 
spiral trajectories are densest around the origin and hence the 
supremum in the definition of path density in (|5]l is achieved 
at the origin. Now let Xi{t) := ctcos(27r(t — i/N)) and 
yi{t) := ctsin(27r(t — i/N)) denote the x and y coordinates 
of the trajectory spi. Then the total length of the trajectories 
located within a disc of radius a centered at the origin is given 
by 

N !L 



Jo 

^iVt'47r2 + 1 



+ — log(i + + l/(47r2)) 



+ o(a2). 



The result follows. 



G. Proof of Theorem \2.8\ 

Before we proceed we need the following lemmas. 

Lemma A.2: For any d > 2 let P denote a uniform set in 
R*^ as defined in (IZSl l. Then for any u < a and any r G W^, 
there is a continuous curve of length 'D^{u, r) + 0(a'*~^) that 
contains the portions of the straight line trajectories from P 
that are located within B!^{r). Here 'D^(u,r) is as defined in 

Proof We prove this claim by induction on d. For d = 2 
the uniform set P comprises a set of parallel hnes and the 
ball B^{r) is a circle for any r S M^. We assume without loss 
of generality that the lines are parallel to the x-axis. In this 
case it is easy to visualize a non-self-intersecting continuous 
curve that passes through the portions of these lines within 
_B2(r). A possible curve is one that starts with the portion of 
the line with the lowest ordinate and then proceeds to portions 
of lines with higher ordinates sequentially, using the portion 
of the circumference of the circle between adjacent portions 
as a connection. Hence the total length of such a curve would 
be no more than the 'D^{u, r) + 2ttu. Hence the statement of 
the lemma is proved for d = 2. 

Now assume that the statement of the lemma is true for 
d = D — 1. Under this assumption we will now establish 
the lemma for d = D. Let r S and w < a be arbitrary. 
Consider a uniform set of the form dZSl l. Let K, :~ {k G 
Z : 3m G Z^~\t G R such that niD-i = k and Pm{t) G 
B^(r)}. Let pk {p™ : m G Z^-\mD-i = k}. Clearly, 
P = {pk : fc G Z}. Now, for each k E K,, let Ik denote the 
intersection of pk with (r) and let £k denote the total length 
of the line segments in Ik- Clearly, J^keK = T>^{u,r). By 
the definition of pk it is clear that the points in Ik lie on a 
hyperplane Hfe in R^ defined by 



Clearly the intersection "Hfe n B^{r) is congruent to a ball in 
R^~^ of the form B?^^ where u <u. Hence by the induction 
hypothesis, it follows that there is a continuous curve Ck that 
covers Ik with length Jk = tk + 0{a^^^). Since the longest 
distance between any two points in B^{r) is no more than 
2a, any two distinct curves Ck and Cj can be joined into a 
single continuous curve of length Ik + ?j + 0(a). Thus there 
is a single continuous curve of length 



E 

fcG/C 



4 + m - l)0{a) = ^{u, r) + \K\{0{a) + 0(a^-')) 



that covers Ukeich- Since > 3, we have 0{a)+0{a^^'^) = 
0{a^^^). Moreover, by the definition of K, it is clear that 
|/C| = 0(a). Thus the total length of the curve is indeed 
'D^{u,r) + 0{a^~^). Thus the lemma is proved by the 
principle of mathematical induction. ■ 
Lemma A.3: Suppose C R'' is a compact convex set with 
a point of symmetry at the origin. Then ilnfi {^'IZi "miUij 7^ 

for some m G Z''"-'^ \ {0} if and only if i "i^i^i € 

Proof: The if part of the statement follows via the 
convexity and symmetry of fl. Now suppose 



/d-i \ 

s G f^nf7 ^ 



for some to G Z^^^^ \ {0} 



If s := s — X]f=i ^i^i^ it follows by symmetry that 
{s, — s, s, — s} C il. By convexity therefore ^{s — s) — 

Now we proceed to the proof of Theorem l2.8l We first prove 
the sufficiency of dZTl ). Suppose P satisfies (|27] |. Lemma lA. 21 
shows that P satisfies condition (C©. Now let us consider 
condition (CHll. Since the field can be measured at all points 
on the lines in P it follows from the structure of P in ( |25] ) 
that we can sample the field at points of the form 

'^m^Vi + mdeVd, m = (toi, TO2, . . . , m^)^ G Z'' (55) 



where e > is arbitrary. The points in ( 155] ) correspond to 
a lattice with finite density in R'' generated by the vectors 
{wi, W2, . . . , Wd} where Wi = Vi,l < i < d — 1 and Wd = 
evd- It follows from fT9\ that the Fourier transform of the 
sampled field, sampled at points on the lattice (1551 ) is made up 
of repetitions of the Fourier transform on the reciprocal lattice 
of points of the form {X^^Li 'TijWi '■ ^ S '^'^} generated by 
the vectors {ui, ti2, . . . , Ud} chosen to satisfy 



{ui,Wj) = 27r(5.y 



(56) 



nk 



kvD-1 + span{ui, U2, . . . , u_d-2, vd}- 



where 5ij denotes the Kronecker delta function. Clearly, the 
first d — I vectors {ui} defined above are exactly as defined 
in the statement of the theorem. Now since e and hence \\wd\\ 
can be made arbitrarily small, it follows that \\ud\\ can be 
made arbitrarily large. Hence for compact sets il, the spectral 
repetitions in the direction Ud can be ignored and hence perfect 
reconstruction is guaranteed provided spectral repetitions at 
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the points {X]f=i : m E ^} do not overlap. Or /. Proof of Theorem \3.2 

equivalently, 



/d-l 



17 n r2 ^ jTZiM, = 0, meZ'^^^\{0}. (57) 



Thus it follows via Lemma |A.3| that under the assumptions of 
Theorem 12.81 we have P e JVn provided 

d-l 

\ ^ niiUi i for all m G Z'^'^ \ {0}. (58) 
1=1 

Now let us consider the necessary condition. Suppose 3m G 

1'^-^ \ {0} such that \ J^t^i ^ h. Consider the field 

defined by 

d—l d I i\ 

/, (r ) = sin( i ^ 771, (r, 7.,) ) n sine ( ),r 



i=l 



where denotes the unit vector along the i-th prin- 
cipal axis. The Fourier transform of is supported 
within a cubic box of width 2e centered at the points 

{5 X]f=i ~5 X^i^Ti^ "^i^i}- the condition on m it 

follows that there exists some e small enough such that 
/e G Bfi- Now, from the condition on {u j given in the 
statement of the theorem, it follows that fe{pm{t)) = for 
all t G R and for all m G where Pm{t) is as defined in 

(l24l i. Hence cannot be recovered from its samples on Pm{t) 
since it is not distinguishable from the field g G Bn that is 
identically - i.e., g{r) = for all r G M''. Thus P ^ Wo. 



From Theorem 12.41 it easily follows that P*^ G £n- Suppose 
P^ is expressed in the form Pi as defined in (fT3] l. Then the 
vector ui defined in (fT4] i corresponding to Pi = P^ is given 
by 



and hence by the definition of u and by iT% we have P^ G 5o 
for all e > 0. Also from Lemma 12.21 we have 



lim.(P^).m, 



In order to complete the proof we will establish that for 
all P G £n, we have £{P) > Suppose P G £n is 

expressed in the form defined in (fTSI l. Using the notations and 
assumptions of Theorem l2.4l it follows that if P G £a we need 



Q ^ n{s) for all s G M'' 

o 

where il{s) denotes the interior of set il(s). Together with the 
relation dSST l satisfied by the width of a convex set, it follows 
that there should exist two points in Q separated by a distance 
greater than Wi^l). This means that 

N 

II ^fcjuj > for some G {-1,0, 1},1 < i < iV. 



H. Outline of proof of Lemma \2.9\ 

Lemma 12.91 can be proved using the same approach as in 
the proof of Lemma l272l Just as we approximated a circle with 
rectangles in the proof of Lemma l2.2l we will now approximate 
a d-spherical baU with parallelotopes. As before, let Pf (x) 
denote a d-dimensional spherical ball of radius a centered at 
X G M''. For r G K'^ let Pi{r) denote the coefficients in the 
basis expansion r = X]f=i Now for m G U^^^ define 

Srn ■■= {r G Biix) : /3,(r) G [m,,m, + 1), 1 < i < d - 1}. 

Let 

Ma := {m G U^^^ : Pr,i{t) G P^(a;) for some t G M} 

with Pm{t) — ^iZi m-iVi + tvd as defined in ( |24] |. For each 
m G Ma let im{a) denote the length of the line segment 
representing the intersection of the line p,n{-) with S^- We 
then have 



Vol(P,^(x)) 



V0l{Sra)+0ia''-'). 



Now we can approximate 5m with the parallelotope generated 
by the vectors {vi,V2, . . . , Vd-i,imia)vd} so that Vol(5m) w 
|det(G)| In fact it can be shown that 



Vol(5m)= |det(G)|i€„(a)+o(a''-i). 

me Ma 

Now J2meM ^rn{o) = T)^ [a^ x) and hence the result follows. 



E 

meMa 



This means that the path density calculated using Lemma 12.21 
satisfies 

N N 
ZTT — ' ZTT — ' 



■i=l 
N 



> ;^llEfc^u^ll> 



27r ■ 



/. Sketch of proof for Proposition \4.2\ 

The result follows via the exact same steps as those followed 
in proving Proposition 12 . 3 1 b v considering the field 



N 



i=l 



{ui,r ~ Wi) 



gc{r) = cexp(-i(s,r)) ]^sin ^ ^ 

JJ^ sinc( — ), r G 



It is easily verified that the field gc vanishes at points on 
the hyperplanes in P for all values of c. Hence gc cannot 
be uniquely identified from its values on P. Moreover, for e 
small enough the field g^. has a Fourier transform supported on 
the set ri. Thus P does not satisfy condition (]V(T) and hence 
PiNf- 
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K. Proof of Theorem 14.51 

We only prove (|45] | since ( |46] | follows directly from Propo- 
sition 132] 

Let P be the manifold set defined in (|43] l with defined 
in dTTT i. Let j^^, ^^27 • ■ • : ^d-il ^ M"* form an orthonormal set 
of vectors spanning the subspace Hi. Suppose that for each 
i we are given samples on the parallel hyperplanes in Pi at 
points of the form {wi + jAihi + J2k=i '^fc^^fe ■ Ji '^fc ^ 
Let denote the Fourier transform of the following impulse 
stream /] of samples from the i-Xh manifold set Pi in (fT6] i: 

d-1 

jGZ,nGZ''-i k=l 

d-1 

k=l 

where 5{.) represents the Dirac-delta function in d-dimensions. 
We have the following lemma. 

Lemma A.4: Let C M'' be a compact set. If e is small 
enough then the spectrum of the sampled field in ( |59] l can 
be expressed in terms of the Fourier transform F of the field 
as 



'Yexp{\{jui,Wi))F{LL> + jui), ujen. (60) 



where Ui are as defined in the statement of Theorem 14.31 

Proof: As before, let . . . , J C M'' form an 

orthonormal set of vectors spanning the subspace Hi. Let ri = 
Aihi and rk = for 2 < k < d. Thus the sampling points 

form a lattice defined by {rk : I < k < d} shifted by Wi. It 
follows from |fT9l that the sampled spectrum of the impulse 
stream in (|59] l can be expressed as 



cxp(i ^(nfcSfc, w,;))P(a; + E rifeSfe) (61) 



fc=l 



k=l 



where Si satisfy {sj, rk) = 2TrSjk, I < j, k < d. This implies 
that si = Ui. Also, making e arbitrarily small ensures that 
||rfc|j — > for k > 2 whence ||sfc|| — > oo for k > 2. Hence 
since ft is compact it follows that for e small enough only 
terms with Ui = for all i > 2 contribute to the summation 
in (|6T) for cj G O. Replacing 7ii with j leads to ( |60] |. ■ 
We now proceed to the proof of (l45T l. Suppose 



Q ^ n{s), for all s e R''. 

We have to verify the condition (]V(l) holds. Let F denote the 
two-dimensional Fourier transform of any field / bandlimited 
to il. Let oj E n he arbitrary. Consider the following set 



N := {n e Z 



N 



Unen} 



(63) 



where U is a d x N matrix with i-th column given by the 
vector Ui = 2^^. It is easy to see via the convexity of Q that 
N is a lattice-convex subset of Z^. Furthermore, condition ( |62] i 
implies that translates of the unit-cell are not contained 
within N. To each element n of N we associate the value 



v{n) ~ F{Ld + Un). Now the conclusion ( |60l ) of Lemma |A4l 
can be restated as 



exp(i(a;,Wj))P*(w) 

~ cxp(i(a; + jui, Wi))F{uj + jui), uj^fl. 

These equations can be evaluated at the points {to + Un : n ^ 
N,l < i < TV} to obtain consistent linear equations in v{.) 
of the form with gl^ = cxp(i(w + Un,Wi))F^{uj + Un), 
where represents the spectrum of the samples taken on 
points on the manifolds in Pi. Applying Lemma IaTI to the set 
N defined in i63[ we conclude that for each lj E il the values 
of the Fourier transform {F{uj+Un) : ?i e N} can be decoded 
using the values {F^{uj + Un) : n E N, 1 < i < N}. Since 
uj G was arbitrary this proves that F{oj) can be recovered 
for all cj G 51 thus completing the proof of i45[ . 



L. Proof of Theorem \4.4\ 

From Theorem l4.3l it easily follows that P'^ e £q . Suppose 
is expressed in the form Pi as defined in ( 1411 1. Then the 
vector ui appearing in Theorem 14 . 3 1 corre spondin g to Pi = P'^ 
is given by 

ew(n)\-^ . 

— ^ u 

27r J 

and hence by the definition of u and by ( l45T l we have P' e 5^ 
for all e > 0. Also from Lemma 143] we have 

w{n) 



Ul 



1 



limMK(P') 



27r 



Now consider any P G £q; expressed in the form of ( l43] ). 
Using the notations and assumptions of Theorem l4.3l it follows 
that if P G £^ we need 

Q ^ h.{s) for all s e M'* 

o 

where fi(s) denotes the interior of set Q,{s). Following the 
same steps as in the proof of Theorem 13.21 we get (P) > 

W(f2) 
27r ■ 
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